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ABSTRACT 


The  feasibility  of  inducing  controlled  movement  of  a  para¬ 
lyzed  limb  by  means  of  applying  electric  currents  is  discussed. 

This  is  preceded  by  a  functional  description  of  muscle  and  its 
electrical  excitability  and  the  musculo-skeletal  system  of  the  human 
arm.  Experiments  are  described  which  were  performed  to  gain  a 
basic  understanding  of  the  mechanical  response  of  the  lower  arm 
due  to  stimulation  of  the  biceps  muscle-  The  dynamic  model  of  this 
single  degree  of  freedom  system  thus  obtained  is  then  used  to 
synthesize  suitable  feedback  control  schemes.  Two  basically 
different  systems  are  discussed:  the  control  pf  the  angular  posi*^ 
tion,  and  the  control  of  the  angular  velocity  of  the  elbow  joint. 

The  velocity  cpntroller  is  found  to  be  superior  and  is  implemented 
in  hardware.  Three  types  of  velocity  controllers  are  compared: 
Open  loop,  proportional  closed  loop,  and  proportional  plus  integral 
closed  loop.  The  proportional  plus  integral  closed  loop  is  found  to 
be  superior.  The  author  also  indicates  that  a  eloped  loop  torque 
control  may  perform  as  well  or  better  and  also  is  simpler  to 
implement. 
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CHAPTER  I 


THE  CONCEPT  OF  AN  ORTHOSIS 
POWERED  BY  STIMULATED  SKELETAL  MUSCLE 

A.  NEED  FOR  ORTHOTIC  DEVICES 

There  are  several  million  people  who  have  potential  need 
for  orthotic  devices.  Their  disabilities  range  from  the  paralysis 
of  a  single  muscle  to  a  state  of  total  disability  51s  in  quadriplegia. 
The  disorder  may  be  the  result  of  poliomyelitis,  neurological 
diseases  such  as  multiple  sclerosis,  Parkinson's  disease,  and 

/CO) 

cerebral  palsy  or  a  spinal  cord  lesion.  Smiths  '  estimated  that 
there  were  100,000  people  living  in  the  United  States  in  1959  who 
had  spinal  cord  lesions.  Moreover,  the  number  of  these  paralyzed 
people  is  on  the  increase  due  to  advances  in  medical  care,  in¬ 
creased  longevity,  population  increases,  the  consequences  of  high¬ 
speed  travel,  etc. 

1.  Significance  of  Level  of  Spinal  Cord  Lesion 

The  degree  of  such  a  patient*s  loss  of  function  of  course 

depends  upon  the  level  of  his  spinal  cord  lesion.  The  levels  of 

innervation  of  the  critical  muscle  groups  are  fairly  well  docu- 

(21) 

mentedv  }  so  that  patients  with  spinal  cord  lesions  can  be  classi¬ 
fied  according  to  the  level  of  the  lesion  from  purely  objective 
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2 


examination  of  remaining  muscle  function.  Nerves  enter  and  exit 
from  the  spinal  cord  at  precisely  determined  points  along  its  length 
between  adjacent  vertebra.  There  are  seven  cervical,  twelve 
thoracic  and  five  lumbar  vertebra.  Lesions  of  thq  spinal  cord  are 
classified  as  C,  T,  or  L,  followed  by  the  number  of  the  vertebra 
(increasing  from  head  to  tail).  When  a  patient  is  listed  as  having  a 
lesion  at  a  given  level,  this  is  the  lowest  level  in  which  all  muscu¬ 
lature  remains  functional.  The  major  muscle  groups  added  at  each 

(40) 

critical  level  are  summarized  below.  The  table  is  representa¬ 

tive  of  the  "average”  patient  with  sharp  lesions  at  one  particular 
level.  An  actual  lesion,  however,  may  span  several  levels  and 
produce  bizarre  combinations  of  lost  muscle  function. 


Deltoid,  Biceps  . . .  .  .  .  .  C~5 

Latissimus,  Serratus,  Pectpralis, 

Radial  Wrist  Extensors  . .  C-6 

Triceps,  Finger  Extensors  and 

Flexors  . . . .  C-7 

Hand  Intrinsics,  Ulnar  side  of  Wrist 

and  Fingers  . . .  T-l 

Upper  Intercostals,  Upper  Back  ......  T-6 

Abdominals,  Thoracic  Extensors  . T-12 

Hip  Flexors,  Quadriceps  . .  L~4 
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2.  Upper  versus  Lower  Motor  Neuron  Lesion 

The  neural  path  from  brain  to  muscle  is  composed  of  two 
neurons  in  series:  the  upper  motor  neuron,  and  the  lower  motor 
neuron.  The  synapse  (connection)  between  the  two  neurons  occurs 
at  the  point  of  exit  from  the  spinal  cord.  For  this  reason  then  a 
spinal  cord  lesion  may  produce  either  an  upper  motor  neuron 
lesion  (UMNL)  or  a  lower  motor  neuron  lesion  (LMNL).  In  the 
case  of  an  UMNL,  the  contractibility  of  the  muscle  is  preserved 
for  years  after  the  accident  but  due  to  disuse  of  the  muscle  it 
may  weaken  considerably.  For  the  muscle  with  a  LMNL,  the 
contractibility  is  also  maintained,  but  for  a  much  shorter  time, 
and  unless  the  muscle  is  artificially  exercised,  it  gradually  wastes 
away  and  becomes  replaced  by  inexterisible  collagenous  tissue. 
There  are  reports^^,  however,  that  indicate  that  the  rate  of 
muscle  atrophy  may  be  decreased  by  a  daily  program  of  electrical 
stimulation. 

A  spinal  cord  lesion  at  a  given  level  then:  will  have  no 
effect  on  muscles  which  receive  their  innervation  from  levels 
above  the  lesion,  produce  UMNL  and  LMNL  at  the  level  of  injury, 
and  produce  UMNL  for  all  muscles  which  receive  their  inner¬ 


vation  below  the  lesion. 
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B .  THE  HUMAN  ARM 

The  human  arm  was  chosen  as  the  object  of  this  study  for 
many  reasons.  Among  them  are  the  experience  with  upper 
extremity  orthotics  acquired  by  the  development  of  the  Case 
Research  Arm  Aid,  and  the  absence  of  large  load  bearing  require¬ 
ments  which  are  important  in  lower  extremity  orthotics.  This 
approach  to  the  orthotic  problem  is  concerned  with  the  electrical 
activation  paralyzed  muscle  to  produce  limb  movements.  It 
therefore  becomes  necessary  to  understand  the  functional  anatomy 
of  the  limb . 


I.  Anatomy  and  Kinesiology 

The  skeletal  system  of  the  arm  consists  of  five  principal 
bones  excluding  the  bones  of  the  hand  (Figure  I-l).  Beginning  at 
the  shoulder  girdle  we  see  that  it  consists  of  four  joints.  They 
are  the: 

Sterno- clavicular  joint  -  this  is  the  only  direct 
skeletal  attachment  of  the  ar rp  to  the  body. 

Motion  at  this  joint  is  not  great,  but  it  allows 
elevation  and  rotation  of  the  clavicle  with 
respect  to  the  sternum. 

Acromio-clavicular  -  this  joint  is  similar  to  the 
above  in  construction  and  positions  the  scapula. 

Gleno-humeral  joint  -  this  is  what  is  usually  thought 
of  as  the  shoulder  joint.  It  is  of  the  ball  and 
socket  type  and  permits  movements  about  three 
^txes  of  rotation. 
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Scapulo-thoracic  joint  -  though  not  a  joint  in  the 
strictest  sense,  there  is  no  ligamentous 
attachment,  the  relative  movement  of  the 
scapula  and  thorax  closes  the  kinematic  loop, 
so  to  speak,  of  the  shoulder  girdle. 

Most  movements  of  the  shoulder  involve  articulation  of  all 
the  above  joints.  For  this  reason,  the  center  of  rotation  of  the 
gleno-humeral  joint  wanders  considerably  during  motion. 

The  elbow  joint  is  the  union  of  three  bones:  the  humerous, 
the  radius  and  the  ulna.  The  hum ero-ulnar  joint  is  purely  mono- 
axial,  allowing  only  flexion  and  extension  of  the  ulna.  Pronation 
and  supination  of  the  forearm  occurs  in  a  distributed  sense  along 
the  radio-ulnar  joint.  That  is,  pronation  and  supination  of  the 
forearm  occurs  when  the  distal  end  of  the  radius  swings  around  the 
ulna. 

At  the  level  of  the  wrist,  the  radius,  ulna  and  carpals 
merge  to  form  the  wrist  joint.  The  metacarpals  pnd  phalanges 
then  form  the  hand. 

Skeletal  muscles  are  attached  at  various  points  among  the 
bones  and  permit  the  central  nervous  system  to  control  the  move¬ 
ment  of  the  limb.  There  are  some  fifty  odd  muscles  which  power 
the  human  arm.  Of  these,  approximately  thirty  can  be  stimulated 
by  means  of  electrodes  placed  on  the  surface  of  the  skin-  The 
points  of  electrode  application  are  called  motor  points.  These 
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locations,  for  the  arm,  are  shown  in  Figures  1-2  and  1-3. 

2.  Functions  of  the  Muscles  of  the  Arm 

Skeletal  muscles  can  be  classified  functionally  as:  prime 
movers,  auxiliary  muscles  and  joint  stabilizers.  The  function  of 
the  prime  mover  muscle  is  to  move  a  joint  in  the  direction  of  its 
pull.  The  auxiliary  muscles  on  the  other  hand  either  assist  the 
prime  movers  in  their  function  or  they  may  be  used  to  deflect  the 
pull  of  the  prime  mover  so  that  the  joint  will  perform  a  more 
oblique  movement.  The  joint  stabilizer  is  used  to  maintain  a  con¬ 
stant  joint  angle  in  the  presence  of  counteracting  torque.  This 
upsetting  torque  may  be  due  to  external  loads  or  it  may  be  due  to 
the  contraction  of  a  muscle  which  spans  more  than  one  joint.  Keep 
in  mind  that  although  each  of  the  three  functions  is  distinct,  a 
single  muscle  may  be  called  upon  to  perform  each  function  at 
different  times. 

The  muscles  of  the  arm  were  investigated  in  terms  of  their 
points  of  attachment  to  the  skeleton.  This  was  done  by  studying 
anatomy  texts  and  also  participating  in  a  cadaver  dissection  at  the 
Western  Reserve  University  Medical  School.  The  author  is 
indebted  to  Dr.  Alvin  Freehafer  and  Dr.  Richard  Briggs  for  this 
opportunity.  The  results  of  this  study  was  the  selection  of  easily 
stimulated  muscles  which  could  perform  the  function  of  prime 
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movers  about  the  major  axes  of  rotation  of  the  arm.  The  results 
are  shown  below. 


SHOULDER  JOINT 


Motion 

Muscles 

Abduction 

middle  deltoid 
supra  spinatus 
upper  trapezius 
serratus  anterior 

Adduction 

latissimus  dorsi 

Extension 

posterior  deltoid 

Flexion 

anterior  deltoid 
upper  trapezius 
serratus  anterior 

Internal  rotation 

pectoralis  major 

External  rotation 

infra  spinatus 

Elevation 

upper  trapezius 
rhomboideus  major 

Depression 

gravity 

ELBOW  JOINT 

Extension 

triceps 

Flexion 

biceps 

WRIST  JOINT 

Abduction 

flexor  carpi  radialis 
extensor  carpi  radialis 
longus 

Adduction 

flexor  carpi  ulnaris 
extensor  carpi  ulnaris 

Extension 

extenspr  carpi  radialis 
brevis 

extensor  carpi  radialis 
longus 
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Extension  (cont. ) 


Flexion 


Pronation 

Supination 


extensor  carpi  radial 
ulnaris 

flexor  carpi  radialis 
flexor  carpi 
ulnaris 

pronator  teres 
pronator  quadratus 

supinator 

biceps 


FINGER  JOINTS 

Extension  _ _  extensor  digitorum 

communis 

Flexion  flexor  digitorum  profundus 

flexor  digitorum  sublimus 


C.  PROPOSED  LIMB  ORTHOSIS  POWERED  BY 
STIMULATED  MUSCLE 


The  above  represents  what  are  thought  to  be  the  necessary 
prime  movers  of  the  main  joints  of  the  arm.  In  a  practical  orthosis 
it  is  believed  that  the  functions  of  the  auxiliary  muscles  and  joint 
stabilizers  could  be  assumed  by  the  mechanical  structure.  Since 
the  arm  will  not  be  required  to  perform  any  strenuous  movements, 
the  assistive  function  of  the  auxiliary  muscles  can  be  neglected. 

The  direction  of  movement  of  a  bi^axial  joint  can  be  controlled  by 


means  of  varying  constraints  placed  along  the  axes  of  rotation. 

Such  a  constraint  can  be  provided  by  a  "joint-lock"  with  an  electri¬ 
cal  actuator.  With  such  a  device,  the  movements  of  the  arm  can 
possibly  be  guided  about  these  individual  axes  in  somewhatthe  same 
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manner  that  the  auxiliary  muscles  perform  this  function.  The 
joint  locks  would  also  serve  the  function  of  joint  stabilizers  when 
they  are  in  the  "lock"  state.  Notice  that  the  arm  may  be  held  in 
any  acceptable  position  with  no  muscle  energy  expenditure  by 
passively  locking  all  the  joints.  This  may  be  quite  useful,  especi¬ 
ally  if  the  patient  attempts  fine  finger  or  wrist  movement  about  an 
extended  terminal  position. 

Orthotic  splints  are  sometimes  powered  by  an  external 
power  source  such  as  an  electric  motor  or  a  compressed  gas 
system.  These  systems  are  generally  heavy  if  the  actuator  is 
actually  mounted  on  the  orthosis  or  complex  if  cables  are  used  to 
transmit  the  force.  An  orthosis  which  utilizes  the  existing  mus¬ 
culature  of  the  paralyzed  limb  as  a  prime  mover  would  have  several 
advantages  over  the  conventional  devices: 

1.  The  forces  would  act  on  the  skeleton  just  as  they 
do  in  the  normal  limb,  thereby  avoiding  any 
subluxation  forces. 

2.  There  are  as  many  actuators  available  as  there 
are  usable  muscles. 

3.  The  bulk  of  the  orthosis  should  be  considerably 
reduced. 

4.  The  contraction  and  relaxation  of  the  muscles 
should  stimulate  blood  and  lymph  circulation 
as  well  as  maintain  muscle  bulk. 


CHAPTER  II 


SKELETAL  MUSCLE 

A.  MORPHOLOGY  OF  MUSCLE 

Before  attempting  to  describe  the  functional  behavior  of 

skeletal  muscle,  it  will  be  useful  to  examine  its  structure.  It  has 

been  found,  by  means  of  the  electron  microscope,  that  a  gross 

muscle  is  composed  of  thousands  of  bundles  (fasciculi)  of  muscle 

fibers;  each  fiber  being  10  to  100  microns  in  diameter  and  having  a 

(15) 

mean  length  of  5  centimeters.  Within  each  of  these  fibers 

there  are  several  hundred  thousand  myofibrils,  0.  5  to  10  microns 
in  diameter,  bathed  in  a  viscous  fluid  which  makes  up  about  50 
percent  of  the  bulk  of  the  muscle.  In  turn,  each  myofibril  has 
about  2500  filaments  of  polymerized  proteins,  called  actin  and 
myosin,  arranged  in  an  orderly  fashion  (Figure  II- l)  every  2  to  3 
microns.  A  fibril  of  5  centimeters  in  length  therefore  contains 
about  20,  000  of  these  series  segments  which  are  called  sarco¬ 
meres.  In  cross  section  the  actin  and  myosin  filaments  are 
arranged  in  a  hexagonal  array. 

When  a  muscle  contracts,  the  actin  and  myosin  filaments 
are  forced  to  slide  past  one  another  in  the  longitudinal  direction, 
thereby  shortening  the  muscle.  This  force  is  thought  to  be 
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a)  REST  LENGTH  OF  SARCOMERE 


b)  CONTRACTING 


c)  FULLY  CONTRACTED 


d)  STRETCHED  TOO  FAR  TO  CONTRACT 

FIG.  I- 1  THE  CONTRACTILE  UNIT  OF  MUSCLE - 
THE  SARCOMERE 
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generated  at  a  series  of  points  in  the  region  of  overlap  between  the 
actin  and  myosin  filaments.  These  "bridges M  between  the  two  fila¬ 
ments  are  believed  to  form  spontaneously  when  an  action  potential 
reaches  the  Z-line.  When  the  muscle  relaxes,  energy  is  supplied 
to  break  the  bond  allowing  the  filaments  to  slicfe  freely  past  each 
other.  Figure  II- 1  also  shows  a  sarcomere  in  relative  states  of 
contraction.  Note  that  only  the  S  filament,  which  joins  adjacent 
actin  filaments  is  stretched  during  lengthening.  During  shortening, 
the  actin  filaments  approach  one  another  and  then  overlap  until  the 
myosin  filaments  butt  up  against  the  Z-line.  Each  sarcomere 
therefore  constitutes  a  dynamic  system  with  spontaneous  force 
input  and  an  output  length  which  is  affected  by  a  viscous  drag  and 
the  external  load. 

This  variability  of  contractile  force  as  a  function  of  sarco¬ 
mere  length  causes  the  gross  muscle  to  have  a  similar  character¬ 
istic.  Figure  II- 2  shows  a  typical  length-tension  diagram  of  a 
skeletal  muscle.  The  total  tension  curve  was  obtained  by  stimula¬ 
ting  the  muscle  at  various  lengths  and  recording  the  tension.  The 
passive  tension  curve  was  obtained  by  stretching  the  passive  muscle 
and  recording  tension.  The  difference  between  these  two  curves 
then  must  be  that  tension  which  is  due  to  the  active  component 
alone.  Just  as  for  the  single  sarcomere  it  can  produce  no  tension 
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LENGTH 


FIG.  31  -  2  LENGTH-TENSION  RELATI ONS 
OF  SKELETAL  MUSCLE 
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at  its  extremes  of  length. 

B.  INNERVATION  OF  MUSCLE 

Every  muscle  fiber  has  a  motor  end  plate  somewhere  within 
the  middle  third  of  its  length.  It  is  at  this  point  that  the  fiber 
receives  its  innervation  from  the  central  nervous  system*  A 
single  motor  neuron  supplies  (innervates)  about  100  muscle  fibers 
by  arborizations  which  terminate  at  the  motor  end  plate.  The 
motor  nerve  and  its  muscle  fibers  constitute  the  "motor  unitn  and 
since  muscle  fibers  follow  the  all  or  none  principle,  the  contrac¬ 
tion  of  the  single  motor  unit  represents  the  quantum  of  muscular 
contraction.  To  perform  a  motor  unit  contraction  the  motor  nerve 
supplies  an  impulsive  change  in  electrical  potential  (action 
potential)  to  the  motor  end  plates.  A  similar  action  potential  is 
generated  in  the  muscle  fiber  and  causes  the  motor  unit  to  twitch. 
Let  us  examine  this  action  in  more  detail.  A  nerve  action  poten¬ 
tial  arrives  at  an  endplate;  a  muscle  action  potential  is  generated 
in  the  vicinity  of  the  endplate  and  travels  down  the  fiber  in  both 
directions  at  a  speed  of  about  5  meters  per  second.  Immediately 
following  the  muscle  action  potential  is  a  "wave  of  contraction" 
which  causes  the  muscle  to  twitch.  This  single  fiber  twitch  is 
actually  manifested  by  the  successive  twitches  of  the  Z0,  000 
individual  sarcomeres  alternately  contracting  and  relaxing. 
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These  individual  contractions  summate  in  the  tendon  to  produce 
the  fiber  twitch.  The  motor  unit  twitch  is  then  the  parallel  addi¬ 
tion  of  single  fiber  twitches.  Gross  muscle  tension  is  produced  by 
the  recruitment  of  many  motor  unifs. 

In  normal  muscular  activities  the  muscle  is  often  called 
upon  to  produce  a  steady  tension.  The  central  nervous  system 
produces  such  a  tension  by  feeding  action  potentials  to  the  motor 
units  at  a  rate  which  is  fast  enough  to  prevent  the  muscle  from 
relaxing  between  twitches.  This  type  of  steady  contraction  is 
called  a  tetanus. 

C.  MECHANICAL  PROPERTIES  OF  MUSCLE 

The  mechanical  properties  of  skeletal  muscle  have  been 
under  investigation  for  many  year s .  Not  until  the  1920's,  however, 
did  adequate  instrumentation  schemes  become  available.  It  was 
during  this  time  that  A.  V.  Hill  began  his  lifelong  work  of  investi¬ 
gating  the  energetics  of  muscle  contraction.  Since  then  countless 
articles  have  been  written  by  him  and  many  others  concerning  the 
mechanical,  chemical  and  thermodynamic  properties  of  skeletal 
muscle.  It  would  be  a  mammoth  undertaking  indeed  to  attempt  to 
summarize  these  and  the  many  other  works  in  the  field.  I  will 
choose  rather  to  develop  a  suitable  model  for  muscle  by  drawing 
upon  the  results  in  the  literature  as  needed. 
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Muscle  is  classically  thought  of  as  n  contractile  element 
with  its  own  particular  properties  and  a  series  elastic  element. 

The  series  elasticity  acts  as  an  energy  storage  device  between  the 
prime  mover  and  the  load  and  has  a  smoothing  effect  on  limb  move¬ 
ments .  This  series  elasticity  has  been  found  by  Hill  to  have  a  very 
high  value  of  stiffness  andfor  this  reason  it  allows  the  contractile 
element  to  be  easily  described  e xp erirpen tally . 

1.  Isotonic  Conditions 

That  is,  if  the  muscle  is  loaded  by  applying  a  constant 
force  (isotonic)  to  its  tendon,  and  it  is  then  stimulated,  the  result¬ 
ing  movement  of  the  tendon  is  assumed  tQ  be  the  internal  shorten¬ 
ing  of  the  muscle  since  the  series  elastic  element  does  not  stretch 
appreciably.  Hill,  in  this  way,  arrived  at  an  experimentally  ob¬ 
tained  description  of  the  mechanical  behavior  of  the  contractile 
element  of  skeletal  muscle.  Figure  II - 3  shows  the  plotted  results 
of  such  an  experiment.  The  force  due  to  the  load  (P)  is  plotted 
versus  the  speed  of  shortening  (v).  As  might  be  expected,  the 
velocity  decreases  as  the  load  increases.  The  relation,  however, 
is  in  general  nonlinear  and  is  given  by  the  equation  of  a  rectangular 
hyperbola: 

(v+b)(P+a)  =  (P^+a)  b  =  constant 


(m.  1) 


VELOCITY  (^f=v) 


-b 

-a 


FIG.  H-  3  RELATION  BETWEEN  LOAD  AND  SPEED 
OF  MUSCLE  SHORTENING 
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The  slope  of  the  curve  at  any  point  has  the  units  of  a  viscosity 
(force /velocity)  and  represents  the  internal  viscosity  of  the  muscle. 
Approximating  the  curve  as  a  straight  line,  an  average  value  of 
muscle  viscosity  (C)  can  be  obtained  and  Equation  II.  1  becomes 


P  =  P'  -  Cv  (II.  2) 

(Refer  to  Figure  II— 3 . ) 

A  simple  model  which  mimics  the  mechanical  behavior  of 

skeletal  muscle  consists  of  a  contractile  element  with  a  parallel 

viscosity  (C)  and  elasticity  (K  ),  and  a  series  elasticity  (K  ) 

(Figure  11-4).  The  internal  and  external  shortenings  of  the  muscle 

are  defined  x.  and  x  .  The  load  velocity  curve  of  Figure  II- 3  can 
i  e 

be  obtained  by  preloading  or  afterloading.  In  Figure  II -4 a  the 
model  has  been  preloaded  by  applying  the  load  to  the  resting  muscle 
then  stimulating  the  stretched  muscle.  Large  loads  will  cause  the 
muscle  to  be  at  a  poor  mechanical  advantage  due  to  its  peculiar 
length-tension  characteristics  of  Figure  11-2,  This  difficulty  can  be 
avoided  by  running  the  experiment  in  the  aftgrloaded  condition  of 
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Figure  II-  Here  the  muscle  at  its  rest  length  lifts  a  weight  from 
a  platform.  Since  the  muscle  delivers  its  maximum  force  at  rest 
length,  larger  loads  can  be  used  and  a  larger  portion  of  the  force- 
velocity  curve  can  be  evaluated.  There  is,  however,  a  slight  delay 
before  movement  occurs  which  increases  with  the  lpad  since  the 
muscle  must  develop  a  tension  equal  to  the  load  before  movement 
begins.  Notice  in  both  curves  that  the  velocity  varies  during  the 
movement.  The  maximum  velocity  is  aqtually  used  in  the  force- 
velocity  curve  since  this  corresponds  to  the  time  at  which  maximum 
muscle  force  is  developed  (neglecting  inertial  force).  The  equation 
of  motion  is  obtained  from  (II.  Z) 


(II,  3) 


2.  Isometric  Conditions 

Under  isometric  (constant  length)  conditions  the  overall 
length  of  the  muscle  remains  unchanged.  However,  the  lengths  of 
the  elastic  and  contractile  elements  change.  The  elastic  element 

The  parallel  elasticity  according  to  Hill  produces  no 
compressive  force. 
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3.  Quick  Release  Experiment 

Quick  release  experiments  have  lpeen  performed  to  evaluate 

the  stiffness  of  the  series  elastic  element  (  Figure  II-6).  The 

muscle  is  first  of  all  stimulated  and  preloaded  and  the  external 

length  (x  )  is  noted.  The  load  is  then  abruptly  released  and  the 
e 

muscle  shortens.  Since  the  contractile  element  cannot  change  its 
length  instantaneously,  the  immediate  change  in  length  is  due  only 
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to  the  shortening  of  the  series  elastic  element-  Then, 


at  t  =  0,  K 

s 


(II.  8) 


4.  Quick  Stretch  Experiment 

If  a  step  change  in  length  is  applied  to  the  resting  muscle 
and  tension  is  recorded,  the  tension  will  initially  rise  to  a  high 
value  and  then  gradually  relax  to  a  lower  level  (Figure  II- 7).  The 
force  equations  are: 


P  =  K  (x  -x.)  (II.  9) 

s  e  l 

dx. 

K  (x  -x.)  =  K  x.  t  C  1  (II.  10) 

s  e  i  pi  dt 


Solving  II.  9  and  II.  10  simultaneously  and  applying  a  step  change 

in  x  yields: 
e 


K  2  x 

P  -  K  x  - 

s  e  K  +K 
P  s 


K  +K 


in.  n) 


Note  that  at  t  =  0,  P  =  K  x  .  Just  a?  in  the  quick  release  experi- 

e 

ment,  the  value  of  stiffness  for  the  series  elasticity  is  available. 

It  should  be  noted,  however,  that  a  quick  stretch  produces  much 
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more  artifact  at  t  =  0  than  does  a  quick  release  due  to  mechanical 
shock  to  the  experimental  setup.  Note  also  that  for  large  t, 


P 


K  K 
P  B 

K  +  K 
p  s 


x 

e 


(II.  12) 


Length-tension  diagrams  of  muscle  are  plotted  using  these  equili¬ 
brium  values.  Since  we  are  dealing  with  linear  models,  the  length 

tension  curve  of  the  model  will  not  have  a  parabolic  shape  as  does 

K  K 

n  s 

muscle  (Figure  II-2)  but  will  be  linear  of  slope  — — • 

p  s 


5.  Twitch  Contraction 

Figure  II-8  shows  the  isometric  twitch  response  of  a  maxi¬ 
mally  stimulated  skeletal  muscle.  Notice  that  first  of  all  there  is 
a  delay  between  the  time  of  stimulus  application  until  the  time  that 
the  tension  begins  to  increase.  This  is  called  the  latent  period  and 
is  usually  between  1  and  10  milliseconds.  The  twitch  curve  itself 
has  a  form  which  resembles  the  impulse  response  of  a  system 
which  is  at  least  second  order.  The  dashed  curve  in  Figure  II-8  is 
the  delayed  impulse  response  of  a  second  order  system  which  has 
been  matched  to  the  experimental  curve.  The  equation  describing 


the  dotted  curve  is: 


TENSION 
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FIG.I-8.  TWITCH  CONTRACTION  OF  MUSCLE  AND 
MUSCLE  MODELS 
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Twitch 


Pt'e 


-25t' 


(II.  13) 


this  is  shown  as  a  dashed  curve  in  Figure  II-8.  Clearly  the  curve 
does  not  describe  the  muscle  twitch.  The  discontinuous  rise  of 
this  curve  af  t  =  0  does  not  match  the  muscle  twitch.  A  much 
closer  fit  can  be  obtained  if  another  spring  and  dashpot  are  added 
to  Qur  inadequate  model  as  shown  in  Figure  II-9.  Since  the  muscle 
twitch  response  was  critically  damped  it  is  necessary  that 


C  =  C'  =  1_ 

K  =  K*  a 

s  s 


(II  *  15) 


and  once  again  that  K1  produces  no  compressive  force-  Also 

P 

it  is  necessary  that  no  loading  occurs  between  the  primed  and 


unprimed  stages,  i.e., 
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FIG. m- 9  IMPROVED  MUSCLE  MODEL  UNDER  VARIOUS  TEST  CONDITIONS 
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x\  »  x.  (11.16) 

1  i 

The  isometric  response  of  the  model  to  an  impulse  of  force  P^  is 
then: 

2  —at  /  \ 

Twitch  of  model  =  K  x.  =  P  a  te  (11.17) 

si  1 


which  is  identical  to  Equation  II.  13. 

6.  E  f  f  e.c  t  o  f  :  T  r:ah  s  du  e  e  r  pn  Isometric  Muscle  Twitch 

Equation  II.  17  is  the  equation  for  a  twitch  contraction  of 
the  model  which  would  be  obtained  using  a  perfect  tension  trans¬ 
ducer.  That  is,  a  transducer  of  infinite  stiffness  and  damping 
and  zero  mass.  Any  realizable  transducer  therefore  would  affect 
the  form  of  the  twitch  contraption.  JLiet  us  first  of  all  observe  the 
effects  of  a  stiff  transducer.  Figure  IIr9a  shows  the  transducer 
mounted  in  series  with  the  model.  The  stiffness  of  the  transducer 
(kt)  in  series  with  Kg  produces  a  resultant  series  elasticity  of 


K' 


K  K 
P  T 

K  +Km 
s  T 


(II.  18) 


which  changes  the  form  pf  twitch  contraction  to 
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The  force  which  the  muscle  produces  is  from  Equation'll.  17 

(assuming  x  is  small,  i.  e.  ,  K  large) 

6 


p 


(11.23) 


Equating  II.  22  and  II.  23, 

x  2 

=  - - - - - - -y  (II.  24) 

P1  (MTs2+CTs+KT)(s+a) 

p,KTa2  PKt 

P  =  - 5- - —  =  - - -  (11.25) 

(Mts  +CTs+KT)(s+a)  (Mts  +Cts+Kt) 


The  tension  as  measured  by  the  transducer  (P^,)  then  is  affected 
by  an  attenuation  and  a  second  order  lag.  We  will  have  recourse  to 
use  this  expression  in  the  experimental  section  of  the  thesis. 

7.  Effect  ;o  f  :  Transducer  on  Isotonic  Muscle  Twitch 

Let  us  first  of  all  calculate  the  isotonic  response  of  the 
model  to  a  constant  inertialess,  undamped  load  (Figure  II-9c) . 

Again  we  make  the  assumption  that  the  series  elasticity  has  a  large 
stiffness  such  that  x  =  x. .  The  equations  of  motion  for  the  muscle 


model  are: 
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P,  =  C'i'.  +  K'  (x'.-x.) 

1  1  S  1  1 


(11.26) 


P  =  K'  (x>x‘.)  +  C'x. 

S  1  1  X 


(11.27) 


Solving  these  two  equations  simultaneously  and  making  use  of  the 


s"  notation: 


K 1  P  -P(C's+K'  ) 
si  s 


1  (Cs+K'  )(C's+K'  )-K' 
s  s  s 


(II.  28) 


Many  so-called  isotonic  experiments  are  performed  using  weights 
as  the  load.  What  are  the  effects  of  this  added  inertia  on  the 
isotonic  contraction  of  the  muscle  model  ?  (Figure  IIi- 9d).  The 
equations  of  motion  are: 


P  =  C'x'  +  K'  (x'.-x.) 
1  1  s  1  1 


(II.  29) 


O  =  K'  (x.-x'O  +  Cx.  +  Mx. 

s  1  1  1  1 


(II.  30) 


Solving  these  equations  simultaneously  yields: 


P.K' 

1  s 


1  s[MC's2+(MK*  +CC')s+K(C+C')] 


(11.31) 
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Equation'll.  31  has  a  much  different  form  from  Equation'll.  28  which 
represents  a  true  isotonic  contraction.  It  should  be  mentioned, 
however,  that  muscles  in  the  body  usually  act  against  inertial 
loads,  especially  the  proximal  joints  of  the  extremities.  The 
results  of  the  experiments  with  inertial  load  should  therefore  have 
more  relevance  to  actual  muscle  function  than  do  the  isotonic 
experiments.  It  should  also  be  remembered  that  although  the 
muscle  force  P^  is  shown  as  a  constant  in  the  equations,  it  is  a 
function  of  the  number  of  motor  units  active  and  the  length  of  the 
muscle. 

In  summary  then,  a  viscoelastic  muscle  model,  including 
the  dynamics  of  contraction  is  given  in  Figure  11^*9.  The  twitch 
contraction  of  this  model  due  to  a  single  action  potential  (impulse) 
is 

twitch  =  P^a^te  *  (II.  17) 


where 


C  _  C' 

K  “  K! 
s  s 


During  the  steady  state  contractiop  of  the  muscle  model  to  a  rapid 
stimulus  rate  (tetanus),  the  response  may  be  characterized  by: 
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isotonic 

isometric 

If  the  transient  response  of  the  model  to  a  tetanic  stimulus  is  also 
desired  even  in  the  presence  of  transducer  effects  then  the 
dynamic  equations  of  motion  of  the  particular  configuration  can  be 
written  and  solved. 

D.  THE  SERVOMECHANISM  OF  MUSCLE  CONTRACTION 

We  have  seen  that  each  muscle  fiber  behaves  much  like  an 
nth  order  system  with  a  delay  time.  These  muscle  fibers  are 
divided  functionally  into  two  classes:  The  extrafusal  and  intra 
fusal  fibers.  The  extrafusal  fibers,  which  are  much  more  numer¬ 
ous,  act  as  prime  movers;  whereas  the  intrafvisal  fibers,  inter¬ 
spersed  throughout  the  muscle,  act  as  regulators  of  gross  muscle 
length.  Groups  of  these  intrafusal  fibers  make  up  the  muscle 
spindles.  These  usually  run  the  whole  length  of  the  muscle  and 
act  as  length  sensors.  It  has  been  shown  experimentally  that  the 
muscle  spindle  emits  a  string  of  action  potentials  to  the  central 
nervous  system,  the  rate  of  which  is  proportional  to  the  change 
of  muscle  length.  In  addition  to  these  afferent  signals,  it  also 
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receives  an  efferent  signal  from  the  central  nervous  system  which 
biases  the  operating  point  of  the  transducer.  Figure  II-  10  shows, 
in  simplified  form,  the  neural  connections  to  an  extrafusal  and 
intrafusal  fiber.  A  command  signal  is  generated  in  the  motor 
cortex  which  is  then  transmitted  down  the  spinal  cord  via  an  upper 
motor  neuron.  According  to  Granit  the  signal  is  then  transmitted 
out  of  the  spinal  cord  on  an  a  or  y  lower  motor  neuron.  In  the 
case  of  7  innervation,  the  signal  travels  to  the  spindle  fiber  and 
causes  it  to  contract.  This  biases  the  spindle  and  is  interpreted 
by  the  transducer  portion  of  the  spindle  as  a  length  change.  This 
error  information  flows  back  to  the  spinal  cord  along  the  groups  II 
(length  change)  and  la  (rate  of  change  length)  afferent  fibers. 

After  synapsing  with  interneurons  the  signals  reach  the  a  motor 
neurons  which  then  cause  the  extrafusal  fibers  to  contract  until 
the  desired  length  is  assumed.  This  feedback  path  is  known  as  the 
reflex  loop  and  is  responsible  for  the  nknee-jerkM  commonly 
observed  in  a  clinical  examination.  That  is,  when  the  tendon  is 
tapped  the  muscle  spindles  are  stretched  briefly  and  vigorously. 

A  burst  of  action  potentials  then  traverses  the  reflex  loop  and 
causes  the  muscle  to  contract  a  short  time  later.  The  y  innerva¬ 
tion  is  thought  to  be  used  by  the  central-nervous  system  in  the 
execution  of  precise  motions.  For  quick,  less  accurate  motions 


MOTOR 
CORTEX 
AREA  4 , 


FIG.  H- 10  NORMALLY  INNERVATED  SKELETAL  MUSCLE 
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The  step  response  of  this  model  also  shows  the  same  type  of 
response. 


f  (t)  =  X  (K  +  K  e"t/T) 
s  ■  m  1  2 


(II.  35) 


Stark^  Abased  his  model  for  the  muscle  spindle  on  the  visco¬ 
elastic  properties  of  the  fiber.  The  viscoelastic  model  is  shown 
in  Figure  11-11.  It  consists  of  three  elements  in  series:  the  ten¬ 
don  stiffness ,  the  lumped  stiffness  and  viscosity  of  the  intrafusal 
fibers,  and  the  stiffness  and  viscosity  of  the  nuclear  bag  of  the 
spindle.  If  the  equations  of  motion  for  the  model  are  written,  a 
transfer  function  is  pbtained. 


f 


Ax  G  (t,  s  +  1) 

c  _  mv  1 

xw  T  (t  s  +  1)(t  s+1) 

m  £  j 


(II.  36) 


where  Ax  is  the  change  in  length  of  the  nuclear  bag  which  deter¬ 
mines  the  firing  rate.  If  the  proper  assumptions  are  made,  the 
time  constants  and  gain  are  given  by: 


(11.37) 


FIG.  E-  II  SPINDLE  MODEL  OF  STARK 
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(II.  38) 

(II.  39) 

(11.40) 


This  transfer  function  represents  a  lead-lag  network.  Stark  found 
from  his  frequency  response  experiments  that  the  spindle  also  has 
this  lead-lag  behavior.  The  step  response  of  this  transfer  function 
is  given  by: 


f  =  X  G  {  1  +  - - 

s  mm  T3T2 


-t/ t  ~  -t/r 

Tl)  e  "(T2"Ti^e  1}  (n-41) 


Again  it  is  of  the  same  form  as  equations  II.  33  and  II.  35  except 
that  it  has  two  decaying  exponentials. 

Another  element  which  is  important  in  the  control  of  mus¬ 
cular  contraction  is  the  Golgi  tendon  of  the  muscle.  If  the  tension 
in  the  muscle  exceeds  a  preset  level,  the  tendon  organ  informs  the 
central  nervous  system.  It  has  no  efferent  input. 

There  are  other  functional  elements  involved  in  the  control 
of  muscle  such  as  joint  receptors,  pressure  receptors,  etc. 


44 


Unfortunately,  not  a,s  much  is  known  about  these  receptors  or 
exactly  how  they  influence  muscular  contraction  at  the  present 
time.  We  will  therefore  limit  the  discussion  to  the  above- 
mentioned  elements, 

2.  System  ; Operation 

The  simplest  skeletal  joint  is  the  mono-axial  joint  bridge 
by  two  antagonistic  muscles.  A  block  diagram  showing  the 
agonistic  action  of  one  of  these  muscles  is  shown  in  Figure  II- 12. 
The  central  nervous  system  innervation  enters  the  motor  neuron 
pool  and  synapses  with  a  y  or  a  a  motor  neuron.  For  y  inner¬ 
vation  the  train  of  action  potentials  is  transformed  to  a  length 
signal  by  the  D/A  converter.  This  reference  input  (  -Xr)  then  is 
gummed  with  the  muscle  length: 


X  =  X  -  X 
e  m  r 


(II.  43) 


MUSCLE  SPINDLE 
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In  the  event  that  the  required  muscle  force  (F^)  becomes 
large  enough  to  cause  damage  to  the  tendon,  the  Golgi  tendon  organ 
fires  and,  after  a  delay  (r  )  has  an  inhibitory  effect  on  the  a 
innervation  of  the  agonist  and  a  facilitory  effect  on  the  antagonist. 
This  protective  role  is  the  only  known  function  of  the  Golgi  tendon 
organ. 

The  above  paragraphs  are  by  no  means  the  last  word  in  this 
highly  complex  business  of  understanding  the  neurologic  control  of 
muscle.  This  area  is  currently  under  the  investigation  of  many 
researchers,  and  new  developments  and  new  theories  are  forming 
at  a  rapid  pace.  The  purpose  of  this  section  is  only  to  review 
briefly  the  general  nature  of  muscular  control.  Some  of  the 
material  will  have  relevance  in  Chapter  III  where  the  behavior  of 
electrically  stimulated  muscle  will  be  discussed. 


CHAPTER  III 


THE  ELECTRICAL  EXCITABILITY 
OF  MUSCLE 

A.  THE  STRENGTH- DURATION  CURVE 
1 .  Mathematical  Modeling 

The  strength-duration  (S-D)  curve  is  a  graphical  representa¬ 
tion  of  the  excitability  of  a  single  axon,  whole  nerve,  or  muscle  in 
response  to  a  pulsed  DC  electrical  stimulus.  The  strength- 
frequency  (S-F)  cui^ve,  similarly,  is  a  measure  of  the  excitability 
of  nervous  tissue  to  an  AC  electrical  stimulus.  We  will  be  con¬ 
cerned  here  only  with  the  S-D  curve.  Figure  III—  1  shows  the  form 
of  a  typical  S-D  curve  which  we  will  assume  is  for  a  single  axon. 

The  ordinate  can  either  be  stimulus  voltage  (V)  or  current  (I) 
depending  upon  which  variable  is  controlled,  and  the  abscissa  is  the 
time  duration  (D)  of  the  pulse.  The  S-D  curve  represents  the 
smallest  values  of  the  parameters  which  will  cause  nerve  excita¬ 
tion  as  determined  by  the  presence  of  an  action  potential  or  the 
slightest  perceptible  muscle  contraction. 

For  pulses  of  very  large  deration  we  notice  that  the  current 
approaches  a  fixed  level  called  the  rheobase.  As  a  relative 
measure  of  excitability,  the  chronaxie  has  been  defined  as  that 
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FIG.mE-l  STRENGTH-DURATION  CURVE  OF 
EXCITATION 


FIG.  HI- 2  STRENGTH  -  DURATION  CURVE 

PLOTTED  ON  COORDINATES  I  VS.  D 
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pulse  duration  which  corresponds  to  a  current  of  twice  rheobase 
as  shown. 


chronaxie  =  Di  (III.  l) 

4  or  V  =  ZxRheobase 

Slightly  different  values  for  chronaxie  are  obtained  depend¬ 
ing  upon  whether  voltage  or  current  stimulation  is  employed.  The 
first  attempt  to  mathematically  model  the  S-D  curve  was  to 
assume  that  the  curve  was  a  hyperbola  In  the  first  quadrant  given 
by  the  equation: 

(I  or  V)  D  =  constant  (III-  2) 


If  we  consider  the  current  stimulus,  this  is  equivalent  to  assuming 
that  the  axon  behaves  as  a  pure  capacitance  and  that  excitation  is 
purely  a  function  of  the  impressed  electrical  charge,  i.  e.  , 


Q  = 


constant 


(III.  3) 


It  is  clear,  however,  that  this  equation  which  may  represent  the 
experimental  curve  within  certain  ranges  does  not  hold  for  large  D 
since  the  experimental  curye  is  asymptotic  to  the  rheobase,  not 
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the  zero  stimulus  axis.  A  better  approximation  therefore  is  the 
equation: 


(I-l  )•■  D  =  constant 


(III.  4) 


where  I  =  rheobase  current. 

R 

This  equation  was  first  proposed  by  Weiss  in  1901  in  the 


form: 


I  =  1R  +  |  (HI.  5) 

Obtaining  the  value  of  MI  M  from  the  experimental  results  is  quite 

straightforward.  The  value  of  "a"  may  also  be  obtained  quite 

easily  from  the  S-D  curve  by  noting  that  when  the  stimulating 

current  is  set  equal  to  21  ,  the  pulse  duration  (chronaxie)  is  given 

R 

by: 

D  =  (III.  6) 

lR 

This  procedure  for  evaluating  the  constants  in  Equation 
(III.  6),  however,  gives  no  indication  of  the  goodness  of  fit  of  this 
equation  to  the  data  points  unless  some  theoretical  points  are  also 
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plotted  and  compared.  A  simpler  way  of  evaluating  the  constants 
and  observing  the  goodness  of  fit  is  to  plot  the  data  on  the  coordin¬ 
ates  I  vs  1/D  (Figure  III- 2).  On  these  coordinates  the  data 
points  should  fall  on  a  straight  line  with  vertical  intercept  "1^" 
and  slope  "a". 

In  some  instances  the  data  points  may  be  adequately  de¬ 
scribed  by  Equation  (III.  6),  but  many  times  the  data  points  do  not 
fall  on  a  straight  line.  A  closer  fit  in  these  instances  can  be 
obtained  by  a  more  geperal  equation  of  the  form: 


(III.  7) 


The  three  parameters  (I  ,  a,  and  b)  may  be  evaluated  by  first  of 

-R 

all  manipulating  the  equation: 


(III.  8) 


log  (I-IR)  =  log  a-b  log  D 


(III.  9) 


This  is  the  equation  of  a  straight  line  on  a  log-log  plot  (Figure 
III-3).  The  value  "b"  is  the  slope  of  the  line. 


53 


FIG.  nr -3.  STRENGTH -DURATION  CURVE  PLOTTED  ON 
LOG -LOG  COORDINATES. 


EXTERNAL 


FIG.  3H- 4.  HODGKIN -HUXLEY  MODEL  OF  SQUID 
AXON  MEMBRANE . 
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m 

x 

m 

y 


(III.  10) 


why  Ay  and  Ax  are  linear  dimensions,  and  and  m^  are 

the  scale  moduli  of  the  log  scales  (i.e.  ,  linear  dimension  of  one 
cycle). 

The  parameter  ,fa"  is  evaluated  by  noting  that  at  D  =  1,  log  D  =  0, 
and 


iog  (i-ir)  =  iog a  (m.  ii) 

I  -  I  =  a  (III.  12) 

1\ 


A  better  fit  can  be  obtained  to  the  experimental  data  using  this 
equation  since  it  is  more  general.  What  this  approach  lacks, 
however,  is  any  physiological  basis.  It  is  a  pure  curve  fitting 
technique  with  some  vague  relation  to  electric  charge. 

A  different  approach  to  describing  the  S-D  curve  was 

(50) 

attempted  first  of  all  by  Blair  and  later  refined  by  Rashevsky. 

In  this  approach  the  S-D  curve  was  synthesized  by  considerations 
of  ion  migrations.  They  assumed  that  there  was  some  leakage  of 
ions  from  the  area  of  stimulation  as  y/ell  as  ion  buildup. 
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Continuing  with  our  electrical  analogy  this  is  equivalent  to  a 
parallel  combination  of  a  resistance  and  capacitance.  We  can  de¬ 
scribe  the  physiological  basis  for  this  assumed  model  if  we  ex¬ 
amine  the  single  axon.  The  functional  model  for  the  squid  axon 
developed  by  Hodgkin  and  Huxley is  shown  in  Figure  III-4. 

The  model  is  characterized  by  a  transmembrane  capacitance  (C) 
in  parallel  with  variable  conductances  for  the  ions  sodium, 
potassium  (g  ,  g  )  and  a  leakage  conductance  (g  ).  In  addition 
the  conductances  are  biased  as  indicated  by  the  sources  (V.^,  V^, 

V  ).  When  the  axon  is  stimulated,  several  electrochemical  events 
L 

transpire  which  result  in  the  production  of  an  action  potential. 

That  is,  when  the  transmembrane  voltage  exceeds  a  threshold 
value,  the  membrane  depolarizes,  thereby  raising  the  conductances 
of  the  sodium  and  potassium  ions.  There  is  then  a  rapid  influx  of 
sodium  ions  followed  by  an  efflux  of  potassium  ions  which  pro- 
duces  an  action  potential.  The  action  potential  then  propagates 
down  the  axon  as  ions  are  exchanged  along  its  total  length. 

An  action  potential  will  be  generated  if  the  transmembrane 
voltage  Y  reaches  a  threshold  value  V  .  If  a  step  change  in 
current  is  applied  to  the  model,  the  transmembrane  voltage  will 
change  in  the  manner  shown  in  Figure  III- 5.  For  this  positive 
change  in  current,  stimulation  will  occur  at  the  cathode.  The  time 
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FIG.  El- 5  TRANSMEMBRANE  VOLTAGE  DUE  TO 

CURRENT  STIMULUS 


STIMULUS  DURATION  (D) 

FIG.  EL-  6  STRENGTH  DURATION  CURVE  PLOTTED 
ON  SEMI-LOG  COORDINATES 
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at  which  this  occurs  (D)  is  the  shortest  time  that  the  stimulus  must 

be  on  for  an  action  potential  to  be  generated.  These  values  of  I 

and  D  therefore  are  a  point  on  the  S-D  curve  for  the  model.  We 

now  have  all  the  information  we  need  to  synthesize  the  S-D  curve  of 

the  model.  To  find  the  rheobase  current  (I  )  we  note  that  for  very 

JK 

large  D  the  requirement  for  stimulation  is 


!RR  -  V  =  VT 


(III.  13) 


V  +V 
T 


(III.  14) 


where 


gNaX  gK 


the  remaining  portion  of  the  curve  is  given  by: 


VT  =  IR(l-e'D/RC)  ,  V 


(III.  15) 


"  V  +V 

T/l  -D/RC,  T  _  T 

I(1-  >  =  *  ht 


(III.  16) 


JE  '  i-e'D/EC 


(III.  17) 
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It  would  be  rather  difficult  to  fit  experimental  data  points  to  con¬ 
form  to  this  requation  in  its  present  form. 

Rearranging: 

'r  ,  -D/RC 

—  =  1-e 

I_I 

R  -D/RC 

T“  =  6 

I  D/RC 

log(-i-)  *  D/RC 
WR 

If  the  data  points  are  now  plotted  on  a  semi-logarithmic 
plot  (Figvire  III-6),  they  should  lie  on  a  straight  line  with  a  slope 
1/RC  =  1/t  where  r  is  the  time  constant  of  the  membrane. 

This  equation  for  the  S-D  curve  fits  the  experimental  data  as  well 
as,  or  better  than,  the  previously  derived  expressions  and  is  to  be 
preferred  since  it  has  some  physiological  significance. 


(III.  20) 

(III.  21) 


(III.  19) 


(III.  18) 
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1*  E  f  f  e  c  t  6 f '  U  sing  Surface  Electrodes 

The  discussion  thus  far  has  been  limited  to  transmembrane 

stimulation  of  a  single  axon.  This  thesis,  however,  is  concerned 

with  the  gross  mechanical  response  of  muscle  due  to  an  electrical 

stimulus  applied  to  the  surface  of  the  skin.  Exactly  how  such  a 

stimulus  causes  membrane  depolarization  is  unknown.  Lane  and 

(35) 

Zebo  have  investigated  the  electric  field  distribution  in  a 
stimulated  cat  limb  in  an  effort  to  determine  if  this  quality  of  the 
stimulus  is  responsible  for  the  induced  contraction.  Among  their 
findings  is  the  fact  that  most  of  the  voltage  drop  occurs  within  the 
skin  just  beneath  the  stimulating  Mcathodeu  electrode.  Further¬ 
more,  this  cathodal  field  distribution  is  relatively  unaffected  by 
variations  in  the  placement  of  a  larger  "anode11  electrode.  They 
are  currently  investigating  the  theory  that  the  gradient  of  the 
electric  field  along  the  length  of  a  fiber  is  responsible  for  its 
depolarization.  Precise  knowledge  of  the  complex  electro- 
physiologic  events  which  transpire  within  a  stimulated  limb  would 
indeed  be  of  interest.  Such  knowledge,  although  it  may  be  of  great 
value  in  determining  optimal  stimulation  schemes,  especially  with 
implanted  electrodes,  is  not  completely  necessary  in  order  to 
control  the  contractiqn  of  skeletal  muscle  by  means  of  surface 
electrodes . 
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There  are  certain  critical  questions  concerning  surface 
electrical  stimulation  which  were  answered  by  means  of  simple 
experiments.  It  is  generally  known  that  muscle  fibers  have  a 
much  larger  chronaxie  than  do  nerve  fibers.  Lane  and  Zebo,  by 
injecting  curare,  were  able  to  block  motor  end  plate  transmission 
and  thereby  obtain  strength  duration  curves  for  muscle.  Muscle 
was  found  to  have  a  chronaxie  of  about  2.  0  ms.  whereas  the 
chronaxie  of  nerve  is  about  0.  2  ms.  When  a  strength-duration 
curve  is  plotted  for  a  normal  muscle,  the  chronaxie  value  con¬ 
forms  to  that  of  nerve.  For  this  reason  then  we  can  assume  that 
surface  electrical  stimulation  produces  contraction  primarily  by 
depolarizing  nerve,  although  some  muscle  fibers  are  probably 
also  stimulated  directly. 

As  we  have  seen  in  Chapter  II,  the  motor  nerve  is  com¬ 
posed  of  many  sizes  of  fibers  both  afferent  and  efferent  which  have 
different  functions .  Most  numerous  among  these  fibers  are  the  a 
motor  neurons.  Is  it  not  possible,  however,  that  the  y  fibers 
are  also  innervated  by  the  electrical  stimulus  thereby  activating 
the  reflex  loop  ?  This  is  believed  not  to  be  the  case  for  several 
reasons.  First  of  all,  the  threshold  of  the  larger  a  fiber  is 
much  lower  than  that  of  the  7  fiber.  The  7  fiber  therefore 
requires  a  stronger  stimulus  in  order  to  fire.  Even  if  a  7  fiber 
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were  activated  along  with  several  a  fibers,  what  would  be  the 
consequences  ?  Let  us  first  of  all  consider  a  twitch  contraction,, 
Since  the  muscle  spindle  produces  negligible  tension,  the  peak  of 
the  twitch  will  remain  unchanged  due  to  spindle  contraction,  but 
the  spindle  may  emit  a  burst  of  action  potentials  since  the  y 
innervation  is  interpreted  as  a  command  to  shorten.  This  burst 
would  then  follow  the  reflex  loop  and  cause  a  second  twitch  of 
lesser  magnitude.  An  experiment  under  the  authors  direction  was 
performed  by  J.  Anderson,  master's  student,  on  a  cat  to  deter¬ 
mine  whether  the  reflex  loop  was  acting  during  stimulation.  He 
measured  tension  of  the  gastrocnemius  muscle  of  the  cat  during 
stimulation  of  its  motor  nerve.  A  twitch  contraction  was  recorded 
with  the  nerve  intact.  The  experiment  was  repeated  after  the 
motor  nerve  was  cut  proximal  from  the  stimulating  electrodes. 
There  was  no  observable  difference  between  the  two  records.  The 
conclusion  was  drawn  from  the  above  experiments  that  surface 
stimulation  innervates  primarily  the  a  motor  neurons  of  the 
stimulated  muscle. 

3.  Clinical  Significance  of  Strength-duration  Curve 

A  standard  diagnostic  measure  of  the  physical  condition  of 
a  skeletal  muscle  is  the  strength-duration  curve  as  obtained  by 
surface  electrical  stimulation.  In  practice,  the  clinician  places  a 
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small,  saline- soaked  cathode  over  the  "motor  point'1  of  the  muscle. 
The  motor  point  is  here  defined  as  that  area  on  the  surface  of  the 
skin  which*  when  stimulated*  produces  the  greatest  contraction 
Of  course,  the  location  of  this  point  depends  on  the  position  of  the 
limb.  A  large  reference  "anode"  is  placed  elsewhere  on  the  body 
to  complete  the  electrical  circuit.  The  clinician  then  systemati¬ 
cally  increases  the  stimulus  amplitude  (voltage  or  current)  and 
duration  until  he  observes  the  least  perceptible  muscle  twitch  and 
records  points  on  an  S-D  curve.  From  this  curve  he  can  then  tell 
whether  the  muscle  is  normally  innervated  or  has  an  upper  or  lower 
motor  neuron  lesion.  Also,  if  he  continues  to  record  S-D  curves 
for  a  period  of  time  he  can  tell  if  the  muscle  is  undergoing  the 
"reaction  of  degeneration".  The  criterion  used  for  these  judgments 
is  the  fact  that  muscle  has  a  much  larger  chronaxie  than  nerve. 
Since  nerve  has  a  lower  threshold  to  electrical  stimulation  than 
does  muscle;  when  stimulation  is  applied  to  the  limb  the  intact 
nerve  will  fire  and  cause  the  muscle  to  contract.  The  S-D  curve  in 
this  case  is  that  of  nerve.  If,  on  the  other  hand*  the  muscle  has  a 
lower  motor  neuron  lesion,  the  motor  nerve  cannot  be  excited  and 
the  S-D  curve  obtained  in  this  case  is  that  of  muscle.  During  the 
reaction  of  degeneration  the  S-D  curve  is  observed  to  shift  to  the 
right  and  upwards,  and  many  times  takes  on  bizarre  forms  which 
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indicate  that  it  is  not  wholly  due  to  nerve  or  muscle,  but  a  combina- 
f  38 ) 

tion  of  both. 

B.  THE  ELECTRICAL.  STIMULUS 

1.  Dis. comfort  Due  to  Stimulation 

Many  electrical  waveforms  repetitively  applied  to  a  motor 

nerve  or  muscle  will  cause  the  muscle  to  contract.  An  early  phase 

of  our  research  was  to  determine  which  of  the  many  possible  wave  - 

forms  was  "the  best.  11  There  were  two  criteria  to  be  satisfied; 

patient  comfort,  trauma  was  to  be  kept  to  a  minimum,  and  minimum 

expended  energy  especially  if  a  portable  stimulator  is  ever  to  be 

realized.  The  response  to  the  four  different  waveforms  shown  in 

(72) 

Figure  III- 7  were  investigated  by  Yodovnik,  et  al.  Waveforms 

a  (sine  wave)  and  b  (alternating  pulses)  are  alternating  currents 
whereas  c  (rectangular  pulses)  and  d  (exponential  pulses)  have  a 
single  polarity,  i.e.  ,  direct  currents.  Four  normal  female  and  six 
normal  male  subjects  were  stimulated  by  placing  a  reference  elec¬ 
trode,  size  3  x  5  inches ,  over  the  biceps  and  an  '‘active11  cathode 
electrode,  size  3 /4 x  1  -  1 /4  inches,  over  the  motor  point  of  the 
extensor  digitorum  commuries.  The  hand  of  the  subject  was  placed 
in  a  relaxed  position  on  a  table  and  the  stimulus  was  increased 
until  the  finger  tips  were  lifted  two  inches  above  the  surface  of  the 
table.  Then  the  stimulus  which  caused  this  movement  was 


DIRECT  CURRENT  ALTERNATING  CURRENT 


FIG.  m-7.  WAVEFORMS  OF  STIMULATION. 


FIG.  Ill- 8.  PAIN  BOUNDARY  FOR  PULSED  DC  STIMULI. 
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recorded.  When  the  subject  expressed  a  feeling  of  pain  or  much 
discomfort,  the  experiment  was  stopped  and  the  values  of  the 
stimulus  were  labeled  npain  point.  n 

For  the  AC  stimuli,  the  range  of  frequencies  investigated 
was  20  to  3000  cps.  It  was  found  that  for  frequencies  above  300  cps, 
the  AC  square  pulses  were  comfortable  whereas  for  frequencies 
below  300  cps  the  stimulus  was  painful.  For  the  sinusoidal  stimu¬ 
lus,  the  pain  point  occurred  at  500  cps. 

For  the  DC  stimuli,  fhere  was  no  difference  observed  be¬ 
tween  the  rectangular  and  exponential  pulses.  Unlike  the  AC  stimuli 
which  can  be  characterized  by  two  parameters  (amplitude  and 
frequency),  the  DC  pulses  require  three  parameters  to  describe 
them  (amplitude,  pulse  duration,  and  pulse  interval).  The  dotted 
line  in  Figure  III-8  represents  the  locus  of  pain  points  for  the  DC 
pulses.  The  area  to  the  left  of  the  dotted  line  represents  the  com¬ 
bination  of  acceptable  parameters. 

R  <  25  +  45.  5  log  D  (ill.  22) 

Comparing  the  AC  and  DC  stimuli,  it  seems  that  one  is  no 
more  painful  than  the  other;  however,  fhe  DC  pulses  require  about 
one  percent  of  the  power  required  for  the  AC  stimuli  for  the  same 
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contraction.  It  was,  therefore,  decided  to  investigate  more  com¬ 
pletely  the  DC  rectangular  pulse  type  of  stimulus.  That  is,  we  were 
interested  in  determining  the  effects  of  changing  the  pulsed  DC 
stimulus  parameters  on  the  muscle  contraction. 

2.  Constant'  Current  1  ve  r  sus  Constant  ■ 

Voltage  Stimulation 

There  are  two  classes  of  stimulators  in  general  use;  the  so- 
called  constant  voltage  stimulator  and  the  constant  current  stimu¬ 
lator.  Before  an  intelligent  choice  can  be  made,  the  impedance  of 
the  load  must  be  considered;  the  load  being  two  electrodes  in  con¬ 
tact  with  the  body.  Figure  III- 9a  shows  a  typical  record  of  the 
voltage  between  two  surface  electrodes  in  the  presence  of  a  con¬ 
stant  current  rectangular  pulse  of  10  ma.  and  0.  2  msec,  duration. 
The  voltage  trace  is  characterized  by  an  exponential  rise  and  fall. 

An  electrical  circuit  model  which  has  the  same  response  consists 
of  a  resistor  (R^)  in  series  with  a  parallel  combination  of  a  resis¬ 
tor  (R_)  and  a  capacitor  (C).  This  lumped  parameter  model,  of 
course,  includes  the  interface  impedance  between  electrode  and 
skin.  By  curve  fitting  techniques  the  following  values  were  deter¬ 
mined  from  Figure  III-9, 


R1 = 

200  Q 

R2  = 

2  K 

C  = 

.  05  uf 
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b)  pseudo-constant  voltage 


FIG.  Ill- 9.  ELECTRODE  VOLTAGE  AND  CURRENT 
DUE  TO  RECTANGULAR  PULSE. 
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In  the  above  case  the  maximum  value  of  the  applied  voltage  was 
20  volts.  In  order  for  a  constant  voltage  stimulator  to  apply  a  con¬ 
stant  voltage  pulse  of  20  volts  to  the  same  electrodes  requires  that 
the  initial  current  be: 

I|t  =  0  '  e7  =  m  -  -lainp  <m-23) 

Clearly  then  a  constant  voltage  stimulator  will  not  deliver  a  square 
voltage  pulse  to  the  surface  electrodes  unless  prohibitive  currents 
are  used.  Figure  III~9b  shows  the  voltage  and  current  traces  of  a 
so-called  constant  voltage  stimulator.  This  does  not  mean  to  imply 
that  voltage  stimulation  is  not  as  effective  as  current  stimulation. 
An  advantage  of  constant  current  stimulation  is  simply  that  it  is 
realizable  using  surface  electrodes  whereas  true  constant  voltage 
stimulation  is  notunless  excessively  large  currents  are  supplied. 

A  further  advantage  of  constant  current  stimulation  is  the  fact  that 
changes  in  electrode  impedance  do  not  affect  the  current  through 
the  electrodes.  In  voltage  stimulation,  changing  electrode  impe¬ 
dance,  whether  due  to  electrodes  drying  out  or  shifting,  affects  the 
amount  of  voltage  that  is  actually  applied  to  the  skin.  All  the 
following  experiments  were  therefore  performed  using  a  constant 


current  stimulator. 
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A  current  driver  circuit  was  added  to  the  output  of  a  com¬ 
mercial  voltage  stimulator.  This  unit  was  designed  to  operate  as  a 
source  of  constant  current  rectangular  pulses  whose  frequency  and 
duration  is  determined  by  the  voltage  stimulator.  Provision  was 
also  made  so  that  the  amplitude  of  the  pulse  train  could  be  modula¬ 
ted  by  an  externally  applied  signal.  Figure  III-10  is  a  schematic 
diagram  of  the  current  driver.  For  constant  amplitude  pulses  the 
modulation  input  is  shorted  and  the  current  level  is  set  by  the 
cathode  variable  resistor.  The  voltage  stimulus  input  can  swing 
from  0  to  250  volts,  it  is  capable  of  turning  the  tube  completely  on 
or  off.  When  on,  the  current  is  determined  by  the  cathode  variable 
resistor  and/or  the  bias  supply  and  the  modulation  input.  The 
static  output  characteristics  of  the  current  driver  are  shown  in 
Figure  III- 1 1 .  With  this  apparatus  then  it  is  possible  to  deliver  to 
the  stimulating  electrodes:  a  constant  current  pulse  with  variable 
repetition  rate,  and  a  pulse  train  which  can  be  amplitude  modulated 
by  a  square,  triangle  or  sine  wave  (Figure  III- 12a,  b,  c). 

C.  THE  MOTOR  POINT 

The  motor  point  has  previously  been  described  as  a  point  on 
the  surface  of  the  skin  overlying  a  muscle  which,  when  stimulated, 
produces  a  greater  contraction  than  any  other  point.  This  is 
admittedly  a  vague  definition.  Work  was  therefore  undertaken  to 
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FIG.3E-  II.  OUTPUT  CURVES  OF  CURRENT  DRIVER 
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FIG.  III-12.  AMPLITUDE  MODULATED  CONSTANT  CURRENT 
STIMULI. 
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more  precisely  define  the  characteristics  of  a  motor  point. 

1 .  Experimental  Equipment 

-~-w  An  experimental  apparatus  was  developed  to  investigate  the 
motor  points  of  the  biceps  and  triceps  of  the  arm.  This  equipment, 
shown  in  Figure  III- 13  is  actually  more  versatile  and  is  more  fully 
exploited  in  a  later  section  of  the  thesis.  It  was  designed  to  hold 
the  arm  immobile  in  the  horizontal  plane  at  shoulder  height.  The 
stimulus  is  applied  via  an  electrode  which  is  spring  mounted  on 
horizontal  and  vertical  tracks  allowing  the  stimulus  to  be  precisely 
located  over  the  surfaces  of  the  upper  arm.  The  electrical  path  is 
completed  by  a  neutral  electrode  held  by  the  hand.  The  elbow 
angle  may  be  indexed  in  5°  increments  by  means  of  a  rack  and 
gear.  The  torque  developed  by  the  stimulated  muscle  may  then  be 
recorded  by  means  of  strain  gages  placed  on  the  cantiliver  beam 
holding  the  rack.  If  the  rack  is  disengaged,  the  elbow  is  free  to 
rotate  and  this  angular  movement  may  be  recorded  by  means  of  a 
potentiometer  on  the  end  of  the  vertical  shaft.  The  first  electrodes 
used  on  the  project  were  the  band-aid-snap  type.  This  electrode 
consists  of  a  band-aid,  a  piece  of  aluminum  foil,  and  a  snap  which 
joins  the  two  and  provides  an  electrical  connector  for  application 
of  a  stimulus.  From  experience  it  was  found  that  this  electrode 
could  not  be  used  directly  on  the  skin  for  fear  of  burning  the  skin. 
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That  is,  unless  the  aluminum  foil  is  coated  with  electrode  paste 
and  pressed  dirmly  on  the  skin  maintaining  a  relatively  constant 
impedance  between  the  foil  and  the  skin  over  the  area  of  the  foil, 
large  current  densities  in  isolated  areas  are  possible  which  can 
actually  burn  the  skin.  An  interface  was  provided  between  the  foil 
and  the  skin  by  inserting  a  saline  soaked  gauze  pad  of  the  same  size 
as  the  foil.  This  worked  well  except  that  sometimes  during  move- 
ment,  the  gauze  pad  would  slip  and  cause  the.  foil  to  come  in  con¬ 
tact  with  the  skin  and  produce  a  burn.  Also,  the  band-aid  proved 
to  be  less  than  ideal  for  keeping  a  constant  pressure  between  elec¬ 
trode  and  skin.  A  better  electrode  is  shown  in  Figure  III- 14.  The 
electrical  lead  is  soldered  to  a  copper  disc  and  a  felt  pad  (saturated 
in  saline)  is  placed  between  the  copper  disc  and  the  skin.  Several 
electrodes  of  this  type  were  made  ranging  in  size  from  3 / 8 1  f  dia¬ 
meter  to  1M  diameter  in  increments  of  1  / 8 11 .  As  greater  accuracy 
became  desirable  the  segmented  electrode  shown  in  Figure  III- 15 
was  developed.  This  consists  of  9  electrically  and  mechanically 
isolated  electrodes.  The  electrodes  are  made  of  felt  strips  1/4"  by 
1  1/4"  long  and  fit  into  teflon  holders  against  conducting  strips. 
Individual  clock  springs  apply  approximately  equal  pressure  to 
each  electrode  segment  thereby  maintaining  good  electrode  to  skin 
contact  even  during  muscle  bulge. 


FIG.  Ill— 1 5 .  SEGMENTED  ELECTRODE. 
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2.  Size  of  a  Motor  Point 

Will  moving  an  electrode  away  from  a  motor  point  cause  the 
torque  to  fall  precipitiously  or  will  it  only  decrease  the  torque  a 
slight  amount?  The  following  experiment  was  designed  to  answer 
this  question  at  least  for  biceps  and  triceps  muscles.  A  3/4n  dia¬ 
meter  electrode  was  placed  over  the  muscle  along  a  midline  of  the 

arm  in  increments  of  one  centimeter,  and  a  stimulus  was  applied 

o 

at  each  location.  During  the  experiment,  the  arm  was  held  at  75 
included  angle  and  the  torque  about  the  elbow  joint  was  recorded. 
The  results  of  an  experiment  with  the  biceps  muscle  are  shown  in 
Figure  III-  16.  Here  three  curves  are  shown  representing  three 
different  current  levels.  Notice  that  the  torque  does  not  change 
greatly  as  the  electrode  is  moved  indicating  that  for  the  biceps 
muscle  the  motor  point  is  not  an  easily  defined  point.  The  triceps, 
on  the  other  hand,  exhibits  a  very  localized  motor  point.  The 
lateral  and  medial  heads  of  the  triceps  were  investigated  and  it  was 
found  that  electrode  movement  of  less  than  a  centimeter  caused 
the  torque  to  fall  drastically.  Typical  torque  values  for  stimula¬ 
tion  at  the  motor  points  were:  2  ft-lbs  for  the  medial  head  and 
1 .  5  ft-lbs  for  the  lateral  head  for  a  stimulus  of  15  ma.  0.2  ms 
pulse  duration  and  repetition  rate  pf  50/sec. 


TORQUE  (FT.  LBSj 
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3.  Does  the  Motor  Point  Move  ? 

Since  the  motor  points  of  the  triceps  were  so  well  defined, 
we  were  able  to  determine  whether  it  remained  in  the  same  loca¬ 
tion  or  whether  it  moved  as  the  elbow  was  flexed  and  extended. 
Figure  III-17a,  b  show  the  results  of  experiments  with  each  of  the 
heads.  The  elbow  angle  was  varied  in  fifteen  degree  increments 
while  the  motor  points  were  stimulated  at  each  angle  and  the  loca¬ 
tion  of  the  motor  point  was  recorded.  In  Figure  III-  12a,  we  see 

that  the  motor  point  of  the  medial  head  has  move^  two  centimeters 

o  o 

as  the  elbow  angle  was  varied  between  45  and  90  .  This  excur¬ 
sion  corresponds  roughly  to  the  rOmovement  of  the  muscle  as  a 
function  of  elbow  angle.  That  is,  assuming  that  the  muscle  has  a 
moment  arm  of  2.  5  cm  about  the  elbow  joint,  the  excursion  of  the 
muscle  would  be: 

r0  =  2.  5  X  -j  -  2cm 
4 

With  the  development  of  the  segmented  electrode,  this 
experiment  was  performed  more  accurately  on  the  biceps  muscle 
of  ten  individuals.  This  was  accomplished  by  indexing  the  arm 
every  ten  degrees .  At  each  of  these  twelve  positions  the  electrode 
segments  were  sequentially  energized  by  operating  the  stimulator 
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in  the  "trains"  mode.  Every  second  a  train  of  six  pulses  (0.  2  msec 
duration,  50 /sec)  was  delivered  by  the  stimulator.  A  manual 
switch  was  then  operated  between  trains  to  direct  the  stimulus  to 
the  different  electrode  segments.  Torque,  stimulus  and  elbow  angle 
were  recorded  continuously  at  each  elbow  angle.  A  typical  record 
is  shown  in  Figure  III- 18.  The  electrode  which  produces  the 
largest  twitch  is  said  to  be  over  the  motor  point.  The  maximum 
excursion  of  the  motor  point  for  the  110°  of  rotation  was  found  to  be 
.  625  to  1. 875  with  an  average  value  of  1 . 0  inch.  Again  this  was 
found  to  correspond  roughly  to  the  r  0  movement  of  the  muscle. 

4.  Effect  of  Electrode  Size 

What  is  the  optimal  electrode  size  for  surface  stimulation  ? 
There  were  two  criteria  which  can  be  used  to  determine  an  optimal 
electrode  size.  The  first  was  patient  comfort.  A  smaller  elec¬ 
trode  introduces  a  higher  current  density  to  the  skin  and  causes 
more  discomfort.  On  the  other  hand,  what  is  the  effect  of  electrode 
size  on  the  elbow  torque?  With  the  arm  held  at  75°,  the  biceps  was 
stimulated  with  electrodes  varying  in  size  from  3/8"  to  1"  in  dia¬ 
meter  in  increments  pf  1/8".  The  results  of  a  typical  experiment 
are  shown  in  Figure  III- 19a.  The  same  experiment  was  repeated 

using  the  segmented  electrodes.  Here  a  much  greater  range  in 

2  2 

electrode  areas  was  available:  5/16  in  to  2-13/16  in  as  opposed 
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1/3/67,  J.L. 
BICEPS 


b )  SEGMENTED  ELECTRODE 

FIG.m-19  ELECTRODE  SIZE  VS.  TORQUE 
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to  about  3 / 1 6  in  to  3/4  in  for  the  circular  electrode.  The  result 
of  this  experiment  is  shown  in  Figure  III- 1 9b .  The  torque  reduces 
by  about  sixty  percent  for  a  nine  fold  increase  in  electrode  size. 

Herein  lies  one  of  the  basic  problems  involved  in  surface 
stimulator,  viz.  ,  how  to  locate  the  stimulus  over  the  motor  point 
during  a  movement.  This  problem  will  be  investigated  in  a  later 
section. 

D.  MAGNITUDE  OF  ELBOW  TORQUE  AS  A 
FUNCTION  OF  STIMULUS 

1.  Muscle  Twitch 

When  a  single  supra-threshold  pulse  is  applied  to  a  motor 
point,  the  muscle  contracts.  If  the  muscle  is  not  allowed  to 
shorten  externally  the  tension  developed  is  called  an  isometric 
muscle  twitch.  The  oscilloscope  traces  of  Figure  III-20  show  a 
series  of  twitch  contraction  measured  as  torques  about  the  elbow 
point.  The  biceps  muscle  was  stimulated  with  a  rectangular  pulse 
of  varying  height  which  also  triggered  the  oscilloscope  trace.  To 
arrive  at  some  indication  of  the  relative  linearity  of  the  system, 
several  records  of  muscle  twitch  were  made  for  different  stimulus 
durations  and  amplitude.  Figure  III- 2 1  shows  the  results  of  one 
series  of  such  experiments  where  peak  torque  of  the  muscle  twitch 
is  plotted  versus  stimulus  amplitude  with  pulse  duration  as  the 
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FIG.  HI-20.  TWITCH  CONTRACTIONS  WITH 
INCREASING  CURRENT. 
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FIGHT  -21  MAGNITUDE  OF  TWITCH  VS.  STIMULUS  CURRENT,  PARAMETER  D 
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parameter.  The  curves  all  have  a  threshold  nonlinearity;  however, 
above  this  level  the  curves  may  be  considered  at  least  piecewise 
linear.  A  linearized  approximation  to  these  curves  is  also  shown 
superimposed  in  Figure  jf II —21*  These  straight  lines  are  given  by: 


Peak  torque  =  (.  25D  +  .  03)(l-10)Ft-Lbs  (III.  25) 


2.  Tetanic  Contractions 

If  a  muscle  is  stimulated  by  a  succession  of  closely  spaced 
pulses  it  is  possible  for  the  muscle  to  maintain  a  steady  contraction. 
This  is  called  a  tetanic  contraction.  Figure  III-22  shows  a  series 
of  oscilloscope  traces  of  tetanic  contractions  for  different  stimulus 
frequencies.  The  traces  were  triggered  by  the  stimulus.  Notice 
that  the  magnitude  of  the  torque  reaches  a  maximum  value  at  a 
stimulus  frequency  of  60/sec.  This  qan  be  seen  more  clearly  in 
Figure  III- 2 3  which  is  a  plot  of  tetanic  torque  versus  the  stimulus 
frequency  with  current  as  the  parameter.  For  stimulus  frequencies 
below  20/sec,  the  torque  was  a  sequence  of  twitches  and  no  readings 
were  made.  From  20  to  50/sec  the  average  torque  level  increased 
but  had  some  ripple  of  decreasing  height.  At  about  50Asec  the 
ripple  was  essentially  gone  and  the  torque  reached  its  maximum 


value. 
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FIG.  Ill- 22.  TETANIC  CONTRACTION  WITH 
INCREASING  FREQUENCY. 
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FIG.  III-23.  TETANIC  TORQUE  VS.  STIMULUS  FREQUENCY 
PARAMETER  I. 
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To  see  exactly  how  the  torque  increased  as  a  function  of 
frequency  between  20  and  50/sec,  several  more  data  points  were 
obtained  in  this  region.  Typical  results  are  shown  in  Figure  III-24. 
The  linearized  apprqximations  to  the  curves  are  given  by: 

Torque  =  (4  +  y  )  log  —  (III.  26) 

If  muscle  is  stimulated  at  higher  and  higher  frequencies  a 
point  is  reached  where  the  muscle  cannot  react  to  every  pulse. 

That  is,  every  other  pulse  will  occur  within  the  refractory  period 
of  the  muscle  and  it  will  be  unable  to  contract  in  response  to  this 
pulse.  To  find  out  at  what  frequency  this  occurs,  the  following 
experiment  was  performed  on  the  biceps  muscle.  A  train  of  two 
pulses  (I  =  15  ma,  D  =  0.4  ms)  with  variable  spacing  was  applied  to 
the  muscle.  Oscilloscope  traces  of  the  resulting  twitches  were 
photographed.  For  a  pulse  spacing  of  200  ms  (5/sec)  the  muscle 
produced  two  distinct  twitches  of  equal  height.  As  the  pulse  spacing 
was  decreased  the  two  twitches  merged  into  a  single  twitch  of 
greater  than  twice  the  magnitude  of  a  single  twitch.  (Figure  III-25). 
When  the  pulse  spacing  reached  4  msec,  the  torque  began  to  de¬ 
crease  until  at  2  msec  (500/sec)  the  twitch  was  identical  to  that  due 
to  a  single  pulse.  The  biceps  muscle  therefore  is  unable  to 
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respond  to  a  stimulus  frequency  of  500/sec  and  frequencies  of  this 
magnitude  should  be  avoided.  Also,  since  no  greater  torques  are 
produced  by  frequencies  above  50/sec  (Figure  III-23),  this  fre¬ 
quency  is  to  be  preferred  at  least  from  a  functional  viewpoint. 

Also  at  50/sec  the  pulse  interval,  20  msec,  is  in  the  comfortable 
range  of  stimulation  (Figure  III-8). 

To  see  how  the  individual  twitches  at  50/sec  add  in  a  devel¬ 
oping  tetanus,  the  biceps  was  stimulated  with  first  one  pulse,  then 
by  two  pulses  space  by  20  ms,  then  three,  £tc.  The  contractions 
were  recorded  each  time  and  are  shown  superimposed  in  Figure 
III-26a.  Figure  III-26b  shows  the  results  of  adding  a  sequence  of 
twitches  separated  by  20  ms.  indicating  how  the  muscle  would 
react  if  its  behavior  were  linear.  Aside  from  the  fact  that  the 
tetanic  torque  levels  are  different,  the  contraction  times  are  also 
different.  This  implies  that  the  individual  twitches  add  in  non¬ 
linear  fashion  to  produce  a  tetanic  torque. 

Figure  III- 27  is  a  plot  of  tetanic  torque  (at  50/sec)  versus  I 
with  D  as  the  parameter.  Notice  that  all  of  the  curves  have  a 
parabolic  shape  suggesting  that  the  data  points  would  plot  more 
linearly  on  a  semi-logarithmic  graph.  Figure  JII-28  displays  the 
same  data  on  the  coordinates:  torque  versus  log  I.  These  data 


can  be  represented  by: 
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Torque  =  7.5  log  — —  ft.  lbs.  (ill.  27) 

T 

Figure  III-29  shows  the  same  data  once  more  with  I  plotted 
versus  D  and  torque  as  the  parameter.  The  lower  curve  (T  =  0)  is 
the  familiar  strength -duration  curve  which  represents  the  threshold 
levels  of  the  muscle.  This  curve  has  been  fitted  to  Hill's  equation 
as : 


I 


3.  5 


1 


-6D 

-e 


(III.  28) 


Looking  at  this  figure  in  a  slightly  different  way  we  see  that 
a  given  torque  may  be  achieved  by  setting  either  I  or  D  and  varying 
the  other  parameter.  That  is,  either  I  or  D  may  be  used  as  a  con¬ 
trol  signal  with  equal  effectiveness. 

In  order  for  stimulated  muscle  to  be  useful  it  must  also  be 
capable  of  maintaining  these  torque  levels  for  a  period  of  time.  It 
is  known  that  the  muscle  fatigues  during  stimulation.  For  this 
reason  all  of  the  above  experiments  included  rest  periods 
between  stimuli  of  at  least  30  seconds,  Records  were  obtained  of 
torque  with  respect  to  tinqte  for  a  constant  stimulus.  The  torque 
level  remained  constant  for  about  15  seconds  before  it  began  to 
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diminish.  This,  of  course,  is  for  norrpal  muscle.  For  paralyzed 
muscle  the  fatigue  rate  may  be  greater. 


CHAPTER  IV 


DYNAMIC  MODELING  OF  MUSCLE  CONTRACTION 

A.  ISOMETRIC  CONTRACTION 
1.  Muscle  Twitch  as  an  Impulse  Response 

Figures  III-20  and  III-22  of  the  preceding  section  show,  in 
addition  to  the  levels  of  torque,  developed  by  the  muscle,  the 
transients  involved  in  developing  these  levels.  That  is,  the  torque 
records  are  preceded  by  a  delay  and  exhibit  a  transient  behavior 
which  is  at  least  second  order.  The  delay  may  be  due  to  a  number 
of  causes,  i.  er  ,  conduction  velocity  of  the  nerve  action  potential, 
motor  end  plate  transmission,  propagation  of  the  muscle  action 
potential  and  its  wave  of  contraction,  slackness  in  the  resting 
muscle,  and  transducer  dead  zone.  The  transient  portion  of  the 
records  actually  include  the  dynamics  of  the  lower  arm  as  well  as 
the  transducer.  We  have  seen  in  Chapter  II  that  transducer  effects 
can  be  significant  in  the  dynamic  behavior  of  muscle.  To  deduce 
the  actual  form  of  the  muscle  twitch  from  these  data  requires  the 
determination  of  the  passive  dynamics  of  the  lower  arm  and  the 
torque  transducer.  Once  this  is  obtained  a  deconvolution  pro¬ 
cedure  can  be  used  to  determine  the  actual  muscle  twitch. 
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Figure  IV-1  shows,  in  block  diagram  form,  the  problem  at  hand. 

Transfer  function  of  arm  support  structure  and  torque  transducer. 

Figure  IV-2  shows  the  transient  response  of  the  arm  support 
structure  and  torque  transducer.  This  record  of  the  deflection  of 
the  torque  transducer  (  6 )  was  obtained  by  releasing  the  arm 
support  from  a  deflected  position.  The  response  of  Figure  1^2  is 
similar  to  the  step  response  of  a  second  order  system  with  transfer 
function. 
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The  values  of  £ 

(65) 

Figure  :IY-2V  ' 


and  co 

n 

and  are: 


can  be  evaluated  directly  frpm 


co  =  105  rad/sec 
n 


Transfer  function  of  arm  support  structure,  torque  transducer 
and  passive  arm. 

The  same  test  was  repeated  with  the  arm  of  a  subject  in  the 
arm  support.  The  transient  response  of  this  combination  is  shown 
in  Figure  IV-3  and  has  a  transfer  function  similar  to  equation  1 
with  £  =  .  25  ,  con  =40  rad/sec.  The  difference  between 


03 


TWITCH 


FIG.  IV -1.  THE  DECONVOLUTION  PROBLEM. 
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FIG.  IV -2.  TRANSIENT  RESPONSE  OF  TORQUE  TRANSDUCER 
AND  ARM  SUPPORT  STRUCTURE. 


FIG.  IV- 3.  TRANSIENT  RESPONSE  OF  TORQUE  TRANSDUCER, 
ARM  SUPPORT  STRUCTURE  AND  PASSIVE  ARM. 
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these  values  and  those  above  must  be  solely  due  to  the  passive 
dynamics  of  the  limb. 

The  equation  of  motion  for  either  system  is: 

T(t)  =  J  6‘  +  C  O  +  K  6  (IV.  2) 

where  J  is  inertia,  C  viscous  damping,  and  K  stiffness.  Evalu¬ 
ating  the  constants  for  the  system  without  the  arm;  the  inertia  is 

given  by  the  cumulative  inertias  of  the  parts  of  the  arm  support. 

4 

This  was  found  to  be: 

2 

J  =  .  023  Ft  Lb  sec 

* 

The  stiffness  of  the  torque  transducer  was  also  calculated: 

K  =  285  Ft  Lb/rad 


The  natural  frequency  then  is  given  by 


(IV.  3) 


* 

See  Appendix. 


106 


which  checks  fairly  closely  to  the  yalue  pf  105  rad/sec  found 
experimentally.  The  viscous  damping  coefficient  is  then 


(IV  .  4) 


=  ,  102  Ft  L»b  sec 


When  the  arm  is  fastened  in  the  arm  support  it  will  add  its  inertia, 
some  damping,  but  essentially  no  stiffness  to  the  dynamics.  The 
elasticity  of  the  elbow  joint  was  found  to  be  of  the  order  of  0.  5 
Ft  Lb /rad  (Figure  .  XY- 4) .  The  inertia  of  the  arm  is  easily  calcu¬ 
lated  using  equation  (IV.  3): 


J  =  023  = -^3  023 

-arm  “n  2  (40)2 


=  .  16  023  =  .  137  Ft  Lb  sec2  =  650  lb  in2 


This  value  is  approximately  twice  the  value  given  by  Corell  7  of 
2 

290  lb  in  for  the  average  size  arm.  The  viscous  coefficient  of 
the  arm  has  a  calculated  value: 


ELBOW  TORQUE  (FT.  LBS.) 
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FIG.Er-4  PASSIVE  TORQUE  VS.  ELBOW  ANGLE 
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Since  we  are  interested  in  determining  the  dynamic  behavior 
of  muscle  contraction  independent  of  the  above  effects,  we  must 
account  for  them  in  our  analysis.  Figure  IV- 5  shows  a  muscle 
twitch  as  recorded.  Knowing  the  combined  transfer  function  of  the 
transducer  and  passive  arm  and  the  measured  muscle  twitch,  we 
can  now  perform  the  deconvolution  proposed  in  Figure  IV-  1.  The 
inverse  problem  is  rather  straightforward.  That  is,  if  the  actual 
muscle  twitch  were  knoym,  the  measured  muscle  twitch  could  be 
evaluated  quite  easily  by  convolution  in  the  time  domain.  A  graphi¬ 
cal  procedure  for  performing  the  inverse  of  convolution  (deconvo- 

(51) 

lution)  was  developed  by  Reswick,  v  -  and  is  described  briefly  in 
the  Appendix.  Using  this  procedure,  the  actual  muscle  twitch  was 
determined.  The  solid  curves  in  Figure  IV-6a,  b  show  a  typical 
measured  muscle  twitch  and  the  graphically  determined  actual 
muscle  twitch.  As  a  check  on  the  accuracy  of  this  solution  a 
standard  numerical  convolution  scheme  was  performed  between 
the  actual  muscle  twitch  and  the  impulse  response  of  the  torque 
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FIG.  IV -6.  RESULTS  OF  GRAPHICAL  DECONVOLUTION. 


Ill 


transducer-passive  arm  combination.  The  resulting  curve  is  shown 
superimposed  (dotted)  on  the  measured  muscle  twitch.  The  curves 
match  quite  well  and  the  impulse  response  of  a  second  order  system 
was  then  matched  to  the  actual  muscle  twitch.  This  produced  a 
transfer  function  for  stimulated  muscle  of: 

Ke" ’ °ls 

Muscle  twitch  =  - * - 5 -  (IV.  5) 

(s+20)Z 

K  -,01s 

Muscle  twitch  =  — — — -  ft-lbs  (IV.  6) 

s  +  40s  +  400 

which  is  also  shown  in  Figure  TV -6b.  Allowing  for  the  stimulus 
threshold  nonlinearity, 

K'(l-IT)e",01S 

Muscle  twitch  =  - y -  ft-lbs  (IV .  7) 

s  +  40s  +  400 

for  I  >  IT 

2.  Sinusoidal  Modulation  of  Stimulus,  Frequency  Response 

As  mentioned  in  Section  III,  the  output  of  the  stimulator  can 
be  modulated  by  an  external  signal  generator.  Frequency  response 
tests  were  therefore  performed  to  check  the  dynamic  transfer 
function  of  Equation  (IV.  2).  The  amplitude  of  the  pulse  train 
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(50/sec,  0.4  ms  duration)  was  modulated  by  a  sine  wave  of  variable 
frequency.  Since  the  pulse  duration  is  quite  small,  the  stimulus 
can  be  represented  as: 


l(t)  =  A+B  sin  2 


00 


5  (t-.  02 n)  (IV.  8) 


n  =  0 


From  sampled  data  theory,  '  the  frequency  spectrum  of 
this  signal  is  the  frequency  spectrum  of  the  modulation  signal  re¬ 
peated  along  the  frequency  axis  every  multiple  of  the  pulse  frequency 
(Figure  IV -7). 

Assuming  for  the  moment  that  the  dynamics  of  stimulated 
muscle  are  represented  by  the  twitch  response  of  Equation  (TV .  2), 
we  then  define  the  linearized  muscle  transfer  function  as: 


Torque(s)  _  Ke 


01s 


- .  01s 


s2  +  40$+  400  (s  +20)2 


(IV.  9) 


In  the  frequency  domain  then  the  change  in  torque  is  given 


T(f)  =  G(f)l(f) 


(IV.  10) 


where  f  is  the  frequency  in  cps 
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Since  the  bandwidth  of  G(f)  is  quite  narrow,  only  the  DC 
and  components  of  the  input  will  affect  the  output  torque 
(Figure  IV-7).  In  the  time  domain  then  the  torque  is  given  by: 


that  is  the  frequency  components  of  the  pulse  train  are  filtered  out 
by  the  muscle  and  the  resultant  torque  is  a  function  of  the  modula¬ 
tion  signal  only. 

The  frequency  response  of  the  biceps  muscle  was  then  ob¬ 
tained  in  the  following  manner.  The  modulation  was  set  at: 


l(t)  =  15  +  5  sin(ZTr  ft)  ma 


(IV.  12) 


and  f  was  varied  from  0,  2  to  3.  0  cps.  The  modulation  signal  and 
the  elbow  torque  were  recorded  on  a  two-channel  recorder.  A 
typical  record  is  shown  in  Figure  IV-8.  The  normalized  amplitude 
and  phase  data  from  two  experiments  are  plotted  versus  frequency 
in  Figure  IV-9.  The  data  points  were  best  fitted  by: 

1 

2 

s  +  80s  +  400 


(IV.  13) 
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This  is  similar  to  the  transfer  function  determined  from  the  twitch 
response  (See  Equation  IV.  7)  except  that  there  is  no  delay  time  and 
the  damping  is  greater.  It  should  be  noted,  however,  that  a  delay 
time  would  have  no  effect  on  the  magnitude  curve  and  a  delay  of 
.  01  sec.  would  produce  a  maximum  additional  phase  shift  of  only 
-11°  at  a  frequency  of  3.  0  cps.  One  would  also  expect  that  the 
viscous  damping  in  a  contracted  muscle  is  greater  than  in  a  relaxed 
muscle,  possibly  accounting  for  the  increased  damping  in  the  trans¬ 
fer  function. 

B .  ISOTONIC  CONTRACTION 
1 .  The  Loading  System 

The  dynamic  model  for  isometric  muscle  contraction  will 
now  be  expanded  to  include  the  isotonic  contraction  of  muscle.  A 
constant  Ipad  system  consisting  of  negator  springs  was  designed 
and  built  under  the  author's  direction  by  Keith  Miller,  M.S.  student. 
This  unit  is  shown  mounted  on  the  experimental  apparatus  in 
Figure  IV- 10.  A  negator  spring  resists  with  a  constant  force 
independent  of  length.  With  this  setup,  constant  loads  with  negli¬ 
gible  inertia  may  be  applied  to  the  elbow  joint  from  0.  5  to  7.  5  ft  lbs 


in  increments  of  0.  5  ft  lbs. 
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2.  Motor  Point  Tracking 

As  stated  in  Chapter  III,  the  motor  point  moves  during 
elbow  rotation  much  as  a  point  on  the  muscle  would  probably  move. 
Any  isotonic  experiments  should  be  performed  with  the  stimulus 
always  over  the  motor  point.  Before  each  isotonic  experiment  was 
performed  therefore  the  movement  of  the  motor  point  was  recorded 
using  the  segmented  electrode  (Figure  III- 15). 

It  was  decided  to  stimulate  using  three  electrodes  at  a  time. 
This  gave  a  contact  area  of  .  94  square  inches  which  was  judged  to 
be  quite  acceptable*  The  motor  point  migration  data  was  reduced 
as  shown  in  Figure  IV- 1 1  where  elbow  angle  is  tabulated  versus 
electrode  groupings .  An  X  in  a  compartment  indicates  that  at 
that  elbow  angle  the  motor  point  is  beneath  that  group  of  three 
electrodes.  We  would  now  like  to  locate  the  stimulus  at  these 
points  during  elbow  rotation. 

A  stimulation  commutator  \y as  designed  and  constructed 

(4) 

under  the  author's  direction  by  Charles  Caldwell,  Ph.  D.  student, 
which  would  ensure  that  the  stimulus  was  always  applied  at  the 
motor  point.  This  consists  of  a  printed  circuit  board  upon  which 
are  mounted  twelve  normally  open  magnetic  reed  switches,  spaced 
radially  every  10°  about  the  center  of  rotation  of  the  elbow 
(Figure  IV- 12).  This  unit  is  stationary.  A  bracket  on  the  arm 
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brace  carries  a  permanent  magnet  which  swings  over  the  reed 
switches  during  elbow  rotation.  The  magnet  and  switches  are 
adjusted  so  that  one  and  only  one  switch  is  closed  at  any  time.  The 
twelve  outputs  of  the  switches  are  connected  to  the  input  side  of  a 
small  patchboard  (Figure  IV-  13),  and  nine  outputs  from  the  patch¬ 
board  are  connected  to  the  nine  electrode  segments.  By  placing 
shorting  pins  in  the  appropriate  patchboard  position,  any  set  of 
electrodes  can  be  activated  at  each  angular  position  of  the  arm. 

The  desired  pattern  previously  obtained  as  in  Figure  IV- 11  is  then 
implemented  on  the  patchboard  before  an  isotonic  experiment  is 
performed. 

3.  Muscle  Shortening  Experiment 

With  the  stimulator  and  electrodes  setup  as  above,  the 
ratcheting  mechanism  is  disengaged  allowing  the  elbow  to  move 
freely.  A  constant  load  torque  is  applied  about  the  elbow  joint 
causing  it  to  assume  the  position  of  extreme  extension.  When  a 
train  of  constant  amplitude  pulses  is  then  applied  to  the  biceps  the 
elbow  flexes  as  shown  in  Figure  IV - 14.  That  is,  the  elbow  angle 
assumes  a  constant  velocity  of  flexion  after  a  delay  and  a  lag. 

The  movement,  of  course,  is  terminated  when  the  joint  reaches 
full  flexion.  In  Figure  IV- 1 5  is  plotted  the  load  torque  versus  the 
steady  state  angular  velocity  of  elbow  flexion  with  stimulating 
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current  as  the  parameter.  In  five  of  ten  experiments,  the  data 
points  could  be  represented  fairly  well  by  straight  lines.  Further¬ 
more,  the  average  slope  of  these  lines  was  found  to  be  0.49  ft-lb- 
sec  with  a  range  of  0.44  to  0.  60  ft-lb-sec.  As  mentioned  in 
Chapter'll,  this  slope  is  actually  a  measure  of  the  internal  damping 
of  the  muscle.  Some  of  the  viscous  drag  must  also  be  atrributed  to 
the  joint  capsule  and  fascia. 

4.  Theoretical  Determination  of  Delay  plus  Lag  Time  (tQ) 

Thus  far,  this  section  has  been  concerned  with  the  dynamic 
modeling  of  the  steady  state  behavior  of  isotonically  contracting 
muscle.  Before  a  steady  state  velocity  is  assumed  by  the  limb, 
however,  there  is  a  delay  plus  lag  time.  The  delay  we  assume  to 
be  the  same  as  for  the  isometric  muscle.  The  lag  is  affected  by 
the  lag  due  to  the  isometric  muscle  and  that  due  to  viscous  drag 
and  the  inertia  of  the  limb.  Assuming  the  lower  arm  to  be  a  pure 
inertia: 


T(t) 

T(s) 

G^s) 


J  'f  +  C$ 
s(Js+C)cj)  (s) 

_  <ft(s)  _  l 

T(s)  s(Js+C) 
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s(^L  S  +  l) 


(IV.  14) 
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where  T^.  is  the  torque  produced  by  the  stimulated  muscle.  The 
elbow  angle  therefore  is  given  by: 


,  s  KK'Ie‘T°S  TLK' 

<f>(s)  =  - 7  "  - 

s  (Qs  +  l)(Ps  +  l)  s(Qs  +  l) 


(IV.  20) 


The  response  of  this  transfer  function  to  a  step  input  was  deter¬ 
mined  using  the  MIMIC  language  on  the  Univac  1107  computer. 

The  response  under  various  load  conditions  is  plotted  in  Figure 
IV^17.  Notice  that  in  addition  to  the  decreasing  slope  of  the 
response  due  to  increasing  load,  there  is  also  a  greater  lag  with 
increasing  load.  It  is  this  latter  characteristic  of  the  response 
which  will  be  used  to  evaluate  tjiis  model.  That  is,  there  are 
several  ways  in  which  the  response  of  the  model  can  be  matched  to 
the  obtained  experimental  data.  Those  methods  are  very  time 
consuming  even  when  the  digital  computer  is  used.  A  simple 
graphical  method  was  devised  which  makes  use  of  the  asymptotes 
of  the  responses  to  determine  the  delay  plus  lag  times. 

None  of  the  isotonic  records  were  observed  to  be  under¬ 
damped  (oscillatory),  therefore  a  transfer  function  for  such  re-; 
cords  would  contain  only  real  poles.  That  is,  the  transfer  function 
for  the  system  under  no-load  conditions  could  be  approximated  as: 
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the  response  of  the  system  to  a  step  change  in  stimulus  is 
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then 


<j)^y  =  KA^t  +  ^(-AqTq  +  A ^  -  Aq  5^^)+  transients  (IV.  30) 


n 


=  KA,  +  KA„(t-T_  -  Y  T  )  +  transients 
1  0V  D  *->  n  ' 

n 


(IV.  31) 


Since  <£(t=0)  =  0,  A^  =  0,  therefore 
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This  is  the  equation  of  a  straight  line  with  t-intercept 

T_  +  V  T  .  Therefore,  to  evaluate  the  sum  of  the  delay  plus  lag 
D  n 

n 

times  of  our  muscle  model,  we  can  simply  extend  the  linear  part 
of  the  angle  response  curve  to  the  tTaxis. 

This  is  true,  however,  only  from  the  no-load  records. 
Introducing  a  load  also  affects  the  t-intercept  as  we  have  seen  in 
Figure  IV -17.  Continuing  with  the  straight  line  approximations, 
the  response  of  the  model,  under  load,  is  given  by: 
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This  is  the  equation  of  two  straight  lines  which  when  summed  gives 
the  steady  state  response  of  </>.  Figure  (IV- 18)  shows  the  graphical 
construction  of  this  response.  The  t-intercept  (t^)  of  this  asymp¬ 
tote  can  now  be  found.  This  point  corresponds  to  the  point  where 
(j>  =  0. 


this  function  is  plotted  in  Figure  IV- 19.  This  graph  expresses  the 
expected  value  of  the  t-intercept  for  our  experimental  data  given 
the  load  and  the  torque  that  the  muscle  exerts  (Kl).  Also  plotted  in 
this  figure  are  experimentally  obtained  points  for  three  different 
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stimulus  currents.  Corresponding  values  of  muscle  torque  were 
obtained  from  the  linearized  load-velocity  curves  and  are  also  indi¬ 
cated.  The  data  points  for  I  =  20  and  25  ma  (3.80,  4.25  ft-lbs)  are 
within  15%  of  the  theoretical  4.  0  ft-lb  curve.  The  15  ma  (2.  5  ft-lbs) 
data  points  on  the  other  hand  have  larger  values  of  t^  than  the 
theoretical  curve  has.  No  explanation  of  this  discrepancy  is  given 
although  this  condition  also  prevailed  in  other  experiments. 

5;  Muscle  Lengthening  Experiment 

If  a  muscle  is  stimulated,  allowed  to  shorten,  and  a  heavy 

load  is  applied,  the  muscle  will  lengthen.  The  rate  of  lengthening, 

just  as  in  the  case  of  muscle  shortening,  is  determined  by  the 

(25) 

muscle  viscosity.  A.  V.  Hill'  '  has  shown,  however,  that  the 
viscosity  is  about  twice  as  large  for  the  lengthening  muscle.  A 
series  of  experiments  were  performed  to  see  if  this  could  be 
observed.  Steve  Diamond,  an  undergraduate  student,  performed 
these  experiments  for  a  junior  lab  project  under  the  author's  direc¬ 
tion.  Three  such  experiments  were  conducted.  The  results  of  one 
of  the  experiments  is  shown  in  Figure  IV-20.  The  slopes  of  the 
curves  ranged  from  0.90  to  1.  16  with  an  average  value  of  1 . 0  ft-lb 
seconds.  This  viscous  coefficient  is  about  twice  the  value  found  in 
the  muscle  shortening  experiments  (0.  49  ft-lb  seconds)  thus 
confirming  Hill's  prediction. 
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CHAPTER  V 


CONTROL  SCHEMES  FOR  STIMULATED  MUSCLE 

A.  POSITION  CONTROL  WITH  TWO  ANTAGONISTIC  MUSCLES 

The  first  attempt  to  control  the  movement  of  a  limb  by- 
electrical  stimulation  was  aimed  at  duplicating,  in  a  crude  fashion, 
the  normal  physiologic  control  of  movement.  We  have  seen  in 
Chapter  II  that  normally  innervated  muscle  has  its  own  built-in 
control  scheme  which  has  at  least  position  and  velocity  feedback. 
The  first  external  control  system  was  designed  to  be  a  closed-loop 
position-control  with  position-feedback.  The  elbow  jpint  angle  was 
to  be  controlled  by  alternately  stimulating  two  antagonistic  muscles 
in  response  to  an  error  signal. 

The  principal  flexor  of  the  elbow  joint  is  the  brachialis 
muscle.  Although  this  muscle  possesses  a  motor  point,  the 
resulting  contraction  due  to  an  external  stimulus  is  less  than  the 
contraction  of  the  biceps  muscle  stimulated  with  the  same  current. 
For  this  reason,  the  biceps  was  always  used  as  the  elbow  flexor  in 
our  experiments.  There  is  a  complication,  however,  in  that  the 
action  of  the  biceps  during  flexion  also  supinates  the  pronated  fore¬ 
arm.  This  problem  was  controlled  by  fastening  the  forearm  in  a 
position  which  was  mid-way  between  supination  and  pronation. 
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The  principal  extensor  of  the  elbow  joint  is  the  triceps 
muscle  which  has  three  heads;  the  long  head,  the  lateral  head,  and 
the  medial  head.  There  is  a  motor  point  associated  with  each  of 
the  heads;  therefore,  a  choice  had  to  be  made  among  them.  The 
long  head  is  the  least  powerful  in  elbow  extension  and  was  dis¬ 
counted.  In  general,  the  lateral  head  produced  the  greatest  con¬ 
traction  and  was  used  in  all  the  experiments.  In  the  following  dis¬ 
cussion  whenever  the  term  flexor  is  used  we  will  be  referring  to 
the  biceps  muscle  and  whenever  the  term  extensor  is  used  we  will 
be  referring  to  the  lateral  head  of  the  triceps  muscle. 

1.  The  Controller 

(73) 

A  controller  was  developed'  '  to  modulate  the  amplitude  of 
stimulus  voltage  to  the  flexor  and  extensor  muscles.  The  command 
signal  (input  position)  can  be  generated  by  a  potentiometer  or  any 
suitable  voltage  source.  A  block  diagram  of  the  over-all  system  is 
shown  in  Figure  V-l.  Chronologically  this  system  was  developed 
by  Dr.  L.  Vodovnik  and  the  author  in  1964.  At  this  time  the  stimu¬ 
lus  voltage -input,  torque  output  characteristics  of  muscle  were 
modeled  as 

“V 

K.  e  (V.l) 

4 

* 


Circuit  diagram  of  controller  is  in  Appendix. 


FIG.  V-L.  BLOCK  DIAGRAM  OF  AGONIST -ANTAGONIST  POSITION  CONTROL. 
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which  has  since  been  improved  to 


(V.2) 


(s+a) 


The  two  are  grossly  equivalent,  however,  if  we  consider  the 

second  order  lag  to  produce  a  delay  of  t  =  —  .  That  is,  the  impulse 

a 

response  of  the  second  order  lag  is 


which  reaches  a  maximum  at 


(V.3) 


£  rte-at]  =° 


(V.  4) 


e‘at  -  ate”at  =  0 
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Therefore  the  analysis  should  not  be  greatly  in  error  if  we  let 


T„  =  T.’fr  — 
0  a 


(V.7) 
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Notice  that  the  controller  portion  of  Figure  V-l  contains  a 
discontinuity  at  e1  =  0.  This  has  been  provided  to  ensure  that  the 
stimulus,  when  on,  is  always  above  threshold  (V  ) .  Otherwise, 
an  unwanted  dead  zone  due  to  sub  threshold  stimuli  would  be  in¬ 
cluded.  The  saturation  level  (B)  is  included  to  limit  the  current 
applied  to  the  muscle.  The  musculoskeletal  dynamics  are  char¬ 
acterized  by  the  extremes  of  flexion  and  extension,  and  <^,  an 
inertia  and  a  viscosity.  The  elasticity  about  the  joint,  being  quite 
small  within  </>  and  <f>  (Figure  IV-4)  was  neglected.  A  theoretical 

l  £ 

analysis  of  the  dynamics  of  the  system  may  now  be  performed.  If 
all  the  controller  components  are  adjusted  so  that  the  system 
operates  in  a  linear  mode,  the  open  loop  transfer  function  for  a 
positive  e'  can  be  approximated  as 


GH  (S) 


K1K2K3K4  e‘T°S 
S  (JS  +  Kg) 


(V.8) 


using  the  mean  values  of  the  coefficients  this  becomes 


-0.05s 

GH  ^  =  s(0.  075s  +  l) 


(V.9) 


A  Nyquist  plot  of  this  transfer  function  is  shown  in  Figure  V-2. 
From  control  theory  the  closed  loop  system  will  be  unstable  since 
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the  curve  encircles  the  singular  point  at  - 1 .  A  means  of  compensa¬ 
tion,  therefore,  is  required  for  the  closed  loop  operation  of  the 
system. 

A  twin-T  network  was  built  to  perform  this  function. 

Placed  in  the  closed  loop  it  decreased  the  gain  of  the  system  at  the 
critical  frequency.  From  the  theoretical  curve  of  Figure  V-2  this 
frequency  is  15  rad/sec  or  2.4  cps. 

2.  Experimental  Findings 

The  open  loop  tests  were  conducted  using  a  low  frequency 
sinusoidal  oscillator  (Hewlett  Packard  Model  202A)  as  the  input 
signal.  Continuous  recordings  of  the  input  signal  voltage  and  the 
angle  of  the  elbow  joint  (output)  were  made  on  a  Visicorder  re¬ 
corder.  The  experimental  setup  is  shown  in  Figure  V-3.  During 
an  experiment,  the  frequency  was  varied  in  steps  from  0.  1  cps  to 
2  cps  for  specific  values  of  input  amplitude. 

The  results  of  the  four  such  experiments  conducted  were 
quite  similar.  For  frequencies  below  0.6  cps  the  records  were 
characteristically  as  shown  in  Figure  V-4.  That  is,  the  elbow 
moves  at  a  speed  of  about  2  rad/sec  in  this  case  to  its  extremes  of 
flexion  and  extension  (saturates)  independent  of  the  input  frequency. 
When  the  frequency  approaches  0.6  cps,  however,  the  output  follows 
an  approximately  sinusoidal  oscillation.  Figure  V-5  shows  a 
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FIG.  V-4  INPUT  VOLTAGE  AND  OUTPUT  ELBOW  ANGLE  VS. 
TIME,  f-0.  3cps. 


■■■■■bibb iiiiiiiii ■■■■■■iiai 
■■■■■■■■■■■■■■■■■■■■!■■■■■!■ 
■■■■■■■■■■■■■■■■■■■■■■■■■■■■i 
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 
■■■■■■■■■■■■■■■■■■■■■■■■■■■■I 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 
■■■■■■■■■■■■■■■■■■■■■■■■■■■■a 
aaaar  ^iiaiaiiinv 
■iar<a  laaaaaaaaar*  naaaaaaiBF 
■■aifla^iiaiiaaaaiaiaiaaaaaia'i 
aaar.aai  laiaaaaaiaaaiaaaaaaaai 
■■■•!■■■  iBBBBBBBrjBBBiiBBBBBiBtB 
■B  iBBBtlBBBBBBB'IBBBI'BBBBBBBIiB 
Bl  BBBII'BBBIBBB  IBB  BA  BBBBBBB  IB 
■f  iBBBB!  ■■■■■■■  IBBBB IBBBBBB IB 
B'lBBBBB  BBBBBBVJBBBBBiVBBBBBI  BB 
la  IBBBBBtlBBBBB’IBBBBBITaaBIBr/BB 
■I  BBBBBBUBBBBB IBBBBB1  BBBBB'IBB 
■FtBBBBBBI  BBBBV  BBBBBBB  IBBBB  ABB 
'iBBBBflaaiaair.BBiaaaa  iibbibbi 
I  IB  ■■■■■■.!■■■' IBB  ■■■BBUBBBFiBBB 
BBBBBBBBl'BBB  ■■■■■■■■■  BBB’IBBB 
MBBBIBI1  III  BBBBBBBBB  1Bil*IB 
d^lIBBaiB  IRrL»4BBBBBfli<ar 
HI '■■■■■■!  "  jBBk IBBBBBB'*’  Jill 
BBB  IBBBBBB  JBBB4BBBBIBV  Mill 
■BBW 1IBBBP iaBBBIliaflBB JBIBI 
BBBB.'BBIV  ilBBBBB.^BIBBF  BBBBBB 
■IBBI  ’'If  JllfllBBI  BIBi  BBBBBBB 
■BBBBki _ ^BBBBBB BBBt  .BBBBBBB 

BBBBBB  BBBBBB  BBBBBB*.  .  BBBBIBBB 
■■BBBBBBBBBBaBBBBBBBB BBBBBBB 
BBBBBB ■■■■■■•■■■■■■■■■ ■■■■■■■ 
BBBBBB BBBB BBBBBBIBBBBIBBBBBBB 
■BBBBBB BBBBBB B BBBBBBB B BBBBBBB 
BIBBBBBBBBBBBaBB BBBB II BBBBBB 
BBBBBBBiaBBBBBBBBBaiaiaBaaai 
BBBB BBBBBBB BBBBBBB BBBBBB BBBB 
BaBBBBBBBBBBBBBBBBBBBBBBBBBB 
BBB BBBBBB BBBBBBBBBBaBBBBBBBB 
I BBBBBBBBII BBBBBBB BRaaaaBBBI 
■III ■■BBBIIIBaBBBIBiaBaBBBBB 
BBBIBII B  BBBBBBBBBB  B BBBBBBB BB 


BltlBBBItlifiKIVIVlIlJIIIRIIItlllllVVIl 


BBBBBB BBBBBB BBBBBBBBBBBBBBBBBBBBB BBBBBB 
■BBBBBB BBaBBBBBBBBBBaBBBBBBBBBBBBBBBBBB 
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
BBBBB BBBBBBB BBBBBBBBBBBBBBBBBBBBB BBBBBB 
BBBBBBB BBBBBB BBBBBBBBBBBBBBBBBBBBB BBBBB 
BBBBBBBBBB BBBBBBaaBBaB BBBBBB BBIBB BBBBBB 
BBBBBBBBBB B BBBBBBBBBB BBBBBBBBBB BBBBBBBBi 
BBBBB BBBBBB^ IBBBIBBBBIir  '■■■■■■■■ IBB  BB 
Ik'BBBB BBBBBB  it  BBBBBBBBBB  jb.  BBBBBBBBBB L«*  IB 
B  BBBBBBB  BBP.BB  VBBBBBBBBV  IB.  1BBBI IBBBB  4Bi  B 
B.  IBBBB  BBBB  IBBi  HBBBBBBBGBBt  BBBBBBBBI.BB  t  « 
Bl 'BBBBIBBB IBBI BBBBIBBB IBBB  BBBBBBB  B’iflBM 
Bl BBBBBBBI BBBB IBBBBBB* BBBB. IBBBBBB ■ JaaBl 
BB  BBBBBBBf  IBBBB.  IBBBBBB  V.BBBB11BBBBBBV  IBBB  I] 
BB  IBBBBBB  JBBBBUBBBBBflf  JBBBB1  BBBBBBLBBB 111 
BBI  IBBBBBB  IBBBB  ft  BBBBBB  IBBBBH  IBBHBH’JBBBBB 
■Bl BIBIBI BIIIBiiBBIBilBBiaa IBBIBB 
BBBBBBBBI /IBB BBB  VBBBBV  BBBBBBilBBBBB ■■■■•■ 
BBB  IBaBB  iaiBBBB  lBBliriBlllflli  aaaif  ■■■■;■ 
■BBI1BBBB  BBBBBBBI  BBBBrl  BBBBBB!  BBBBI J  IBBBB  Bl 
■■Bl BBBl BBBBBBBI BBBB  BIB BBIBB IBBB IBBBBBB 
BBBB  BIBF IBBBB II II  IBBf  BBBBBBBBI  'BIB  IBBB  BBfll 
BBBB  IBB liBBBBBBB.  BflNv^flBBflBBI BBf PT1BB«I 

Baaai  aF'  .v  bbbbbbi  iv  ■  -*.  ibbbbbi  ii  lw  bbbbI 

BBBBB"  BBBBBBBBBB. ’*  IBBi  BBBBBBI'  JlilBtll 
Baaaai iBiiivBiRiifc  laaaaiBBiaia/iBiBi  ■  «■] 
aaaiar  BBRaBk'BBiiar  jiiiii  bbbbbb  jbibb  tail 
BBIBV  iaMBBlIBBIff  JBBBBB1  lIBBBf  fllllM.  «B 
IBB BiB BBBBBBl «BI  BBBBBBBI IBBF iBBBBBBI 111 
k.  BBBBIBIBli  ^  .BBBBBBBBB. ’ll  BBBBIBBB*.  ll 
aim.1  BBBBBB  IB  BBlBB  BIBBBB  BBB*. . .  IBBIBBBBB  •■■ 
bbbbbbbibbbiiiibibbbbibbbbbbbbbibbbbbibI 
B BBBBBBBBBBBBBBBB BBBB B  B BBBBIBBB IB BIB IB ll 
BBBBBBBBBBBBB BBBBBB BBBBBBB BBBB BBBBBB BBBl 
BB ■ BBBBIBBBBBBBBBBBBBB  BBBBBBB  BBB  BIBBIial 
BBBB1BBBBBBBB BBB BBB BBBB IBBB BBBl IIBBB llll 

■■■■■bbbbbbi ■■■■■■■■■■■ ■■■■■bbbbbbbbbibI 
■■■■■■■■■BBB ■■■■■■■■■■■■■■■ BBBBBBBBBB BBB 

BBBBBB BBI IBBBBBB «l BIBB BBBB IBBI IB BBBBB I ll 
■BBBBBBI BIBB I BBBBBB BIBB II BBB BBB IBBBBB III 
■■■■■■■■■■■■■■■■■■■ ■■■■■ BBBBB ■■■■aBBBIB| 

■  ■■■■ BBBBBBIBBIIBBBI ■■■■ IBIBIBBIB  fBBBIBl 


FIG.  V-5.  INPUT  VOLTAGE  AND  OUTPUT  ELBOW  ANGLE  VS. 
TIME,  f-0.  8  cps  . 


typical  record  obtained  at  a  frequency  of  0.  8  cps. 

It  was  found  that  for  each  of  the  four  experiments  conducted 
that  the  output  lagged  the  input  by  180°  at  a  frequency  which  was 
invariably  between  0.  8  and  1.  0  cps.  *  Also  the  open  loop  gain  at  this 
frequency  was  greater  than  1  indicating,  as  suspected  from  the 
theoretical  analysis,  that  the  closed  loop  system  will  be  unstable 
unless  some  means  of  compensation  is  provided. 

For  frequencies  of  about  1. 3  cps  the  output  followed  an 
erratic  variation  of  low  amplitide  which  at  times  followed  a  sub- 
harmonic  of  the  input  frequency. 

As  a  further  confirmation  that  compensation  was  necessary, 

the  loop  was  closed  and  the  compensator  was  omitted.  An  input 

signal  was  then  applied  to  the  controller  and  the  arm  began  to 

oscillate  violently  at  about  1  cps.  The  twin-T  network  was  then 

inserted  and  the  oscillations  ceased.  The  compensator  introduced 

a  disadvantage,  however,  in  that  the  dynamic  response  of  the  closed 

loop  system  was  affected.  Because  of  the  attenuation  of  the  filter, 

the  frequency  response  on  either  side  of  0.87  cps  was  also  affected. 

In  the  static  case  the  angular  position  of  the  elbow  joint  was  con- 

+  o 

trollable  to  within  about  -  5  . 

As  a  concluding  experiment,  the  amplifier  of  the  electronic 
bypass^^  was  used  to  provide  an  external  signal  to  the  controller. 
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EMG  pickup  electrodes  were  placed  over  the  extensor  digitorum 
communis  (EDC)  of  the  left  arm;  then  they  were  connected  through 
the  amplifier  to  the  external  (EXT)  input  jack  of  the  controller. 
When  the  wrist  was  relaxed,  the  flexor  received  the  maximum 
stimulus,  whereas  when  the  EDC  was  contracted  maximally  the 
extensor  received  the  maximum  stimulus.  The  subject  was  then 
able  to  control  the  position  of  his  elbow  joint  by  varying  the  con¬ 
traction  of  the  EDC  of  his  left  arm.  No  quantitative  measurements 
were  made  at  this  time  but  the  subject  was  asked  to  point  to  the 
experimenter's  finger  as  the  experimenter  moved  his  finger  ip.  a 
random  fashion  in  the  plane  of  rotation  of  the  arm.  With  very  little 
training  the  subject  was  able  to  perform  this  function  quite  well  for 
speeds  within  the  normal  range  of  movement. 

B .  BIDIRECTIONAL  CONTROL  WITH  A  SINGLE  MUSCLE 
AND  A  FATIGUE-PREVENTING  LOCK 

A  large  class  of  limb  movements  are  performed  against 
gravity.  These  movements  may  be  performed  by  the  action  of  a 
single  muscle  acting  against  the  gravitational  load.  A  consequence 
of  this  scheme  is  that  effort  must  be  expended  to  maintain  a  con¬ 
stant  joint  angle.  Since  the  muscle  is  subject  to  fatigue  it  would  not 
be  practical  to  expect  the  stimulated  muscle  to  maintain  a  constant 
joint  angle  for  any  length  of  time.  A  proposed  solution  to  this 
fatigue  problem  is  to  install  a  joint  lock  and  stimulate  the  muscle 
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only  when  a  change  of  joint  angle  is  desired.  Use  of  a  joint  lock 
requires  that  the  control  system  is  stable.  For  this  reason  the 
first  past  of  this  section  will  be  concerned  with  determining  the 
necessary  parameters  for  stable,  acceptable  control  system  opera¬ 
tion  irrespective  of  the  application  of  a  lock.  The  lock  and  the 
locking  criteria  will  then  be  discussed  followed  by  the  experimental 
evaluation  of  certain  of  the  control  systems. 

1.  Control  Modes 

The  use  of  a  counteracting  load  torque  in  the  control  system 
requires  that  the  muscle  be  operated  about  a  biased  operating 
point.  That  is,  the  muscle  must  produce  a  torque  equal  to  the  load 
torque  to  maintain  a  joint  angle;  and  fluctuations  of  the  muscle 
torque  about  the  operating  point  cause  the  limb  to  change  position. 

If  the  load  is  variable  an  added  complication  arises  in  that  the 
operating  point  must  shift  to  compensate  for  the  varying  load.  Such 
a  varying  load  torque  would  be  produced,  for  example,  when  the 
limb  rotates  about  a  horizontal  axis.  The  moment  arm  of  the  load 
is  then  related  by  the  cosine  of  the  joint  angle.  We  will  here  con¬ 
sider  the  case  where  the  load  is  constant. 

The  first  control  mode  considered  was  the  closed  loop 
position  control  with  a  switching  element  (Figure  V-6).  The 
characteristics  of  electrically  stimulated  muscle  contain 
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nonlinearities,  i.  e.  ,  dead  ?one,  migrating  motor  point,  angle 
dependence  of  muscle  torque,  nonlinear  gain  and  fatigue.  The 
switching  element  gates  the  stimulus  either  on  or  off  and  bypasses 
the  dead  zone  and  nonlinear  gain.  That  is,  the  stimulus  is  either 
zero  or  some  level  above  the  threshold.  This  system  was  modeled 
on  the  Univac  1107  digital  computer  using  the  MIMIC  simulation 
language.  Figure  V-7  shows  the  step  response  of  this  system  for 
two  different  values  of  dead  zone  (E^)  of  the  switching  element. 

The  loop  gain  and  dead  zone  of  the  system  were  varied  in  an 
attempt  to  obtain  acceptable  performance.  In  each  case  the 
response  contained  a  limit  cycle  of  about  1.0  cps.  Increasing  the 
size  of  the  dead  zone  decreases  the  amplitude  of  the  limit  cycle  but 
in  so  doing  the  controller  becomes  insensitive  to  small  changes  in 
the  input  signal.  Also,  the  load  imposes  a  steady  state  error  on  the 
joint  angle.  For  these  three  reasons  then:  a  limit  cycle,  dead 
zone  for  the  input  signal,  and  high  load  sensitivity,  the  closed  loop 
position  control  with  a  switching  element  was  judged  to  be  un¬ 
acceptable.  If  such  a  discontinuous  controller  was  poor,  perhaps 
a  continuous  closed  loop  position  system  would  be  better. 

The  switching  elemenf  was  replaced  by  the  discontinuous 
plus  proportional  element  shown  in  Figure  V-8.  This  is  the  same 
type  of  element  as  used  in  the  controller  of  Section  V  A.  Any  error 
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signal  produces  a  stimulus  which  is  above  the  threshpld  value  of  the 
muscle.  This  system  was  also  modeled  on  the  Univac  1107  and 
the  step  response  is  shown  in  Figure  V-9.  The  loop  gain  (.  084) 
was  calculated  to  give  stable  operation.  The  critical  gain  was 
found  to  be  .  114  at  a  frequency  of  6.  2  rad/sec.  The  response  of 
the  system,  though  stable,  is  also  very  load  sensitive.  The  load 
sensitivity  may  be  calculated, 
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s  ( .  15s  +  l) 


(V.  10) 
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Applying  a  unit  step  change  in  load  and  using  the  final  value 
theorem: 
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The  closed  loop  position  controller  would  therefore  have  a  static 
transfer  function 


6  =6-23.8  T 

c  r  r  Li 


(V.13) 


The  load  sensitivity  is  very  high.  Compensation  schemes,  however, 
could  theoretically  be  applied  which  would  allow  a  higher  loop  gain 
and  therefore  better  regulation. 

The  next  type  of  control  considered  was  a  velocity  control. 
Since  the  limb  is  to  be  controlled  by  a  hpman,  he  can  guide  the 
limb  to  a  desired  position  by  controlling  its  velocity.  The  simplest 
form  of  this  control  is  the  open  loop  velocity  control  as  shown  in 
Figure  V-TO.  This  system,  of  course,  is  stable  but  just  as  the 
position  control,  it  is  quite  load  sensitive.  That  is,  the  greater 
the  load  the  slower  the  limb  moves  for  a  given  input  signal.  This 
error  in  output  velocity  due  to  the  load  is  given  by 


Tt  . 15s+l 
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Applying  a  unit  sfep  change  in  load  and  using  the  final  value 
theorem 
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Closing  the  velocity  loop  (Figure  V-10)  has  the  effect  of 
decreasing  this  load  sensitivity,  i.  e.  , 
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A  unit  step  change  in  load  gives  the  final  value 
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The  critical  loop  gain  for  stability  was  found  to  be  2AK ^  =  3.15  at 
a  frequency  of  20  rad/ sec.  The  error  in  velocity  is  then  equal  to 


6  -  -  .482  rad/sec/ft-lb 

rc 

t=  00 

Even  greater  load  insensitivity  can  be  achieved  by  using  a  propor¬ 
tional  plus  integral  controller  (Figure  V-10),  The  velocity  erro;r 
due  to  the  load  then  is  given  by: 
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The  critical  gains  for  stability  were  calculated.  For  the  propor¬ 
tional  part  a  gain  of  1. 65  was  critical;  for  the  integral  part  it  was 
11.6.  The  frequency  at  the  critical  point  was  15  rad/sec. 

These  three  velocity  control  systems  were  simulated  on 
the  Univac  1107  computer  using  the  MIMIC  simulation  language. 
The  velocity  outputs  of  each  system  to  a  step  input  of  1. 5  rad/sec 
are  shown  superimposed  in  Figure  V-ll.  The  load  sensitivity  of 
the  open  loop  and  closed  loop  proportional  systems  is  quite 
obvious.  Both  the  closed  loop  systems,  however,  are  oscillatory, 
which  may,  at  first  glance,  seem  unacceptable.  Remember, 
however,  that  although  these  are  velocity  controls,  what  we  are 


LOAD  =  10  FT.  LB. 
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FIG.3C-II  RESPONSE  OF  VELOCITY  CONTROLLERS  TO  STEP  INPUT 
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really  trying  tp  control  is  the  position.  The  integral  of  the  velocity 
curves,  the  position  curves,  are  shown  in  Figure  V-12.  Notice 
that  although  there  are  slight  oscillations,  the  curves  closely 
approximate  ramps. 

2.  Joint  Lock 

One  of  the  biggest  problems  encountered  in  this  project 
was  finding  a  small,  lightweight  joint  lock  which  had  the  necessary 
torque  capacity.  The  lock  should  be  capable  of  maintaining  the 
weight  of  the  lower  arm  plus  a  light  load  in  the  hand.  An  average 

(5) 

forearm  and  hand  weighs  3.  7  lbs.  and  has  a  moment  arm  of  6.  7  in. 
Adding  a  2  lb.  weight  in  the  hand  increases  the  load  to  about  4.  2 
ft-lbs.  The  commercially  obtainable  clutches  and  brakes  with  this 
torque  capacity  were  all  about  4  inches  in  diameter  by  about  2 
inches  high  and  weighed  about  10  lbs.  It  was  obvious  that  such  a 
brake  could  not  be  used  on  an  orthosis  so  work  was  begun  to  design 
and  build  one.  Clearly  any  type  of  friction  lock  would  be  too  large. 
A  ratcheting  mechanism  would  support  the  load  and  be  lightweight 
but  it  could  only  by  indexed  at  discrete  angles  which  are  separated 
by  at  least  5  degrees.  This  corresponds  to  positioning  the  finger¬ 
tips  in  increments  of  about  1. 5  inches. 

A  bidirectional  lock  which  would  theoretically  have  infinite 
resolution  was  designed  and  built.  This  lock  shown  in  Figure  V-13 


LOAD  =  1.0  FT.  LB. 
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TIME  (SEC.) 

FIG.  Y- 12  POSITION  RESPONSES  OF  VELOCITY  CONTROLLERS  TO 

STEP  INPUT 
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works  on  the  roller  wedge  principle.  In  the  lock  state  both  pips  are 
forced  into  the  narrowing  annular  space  between  the  rotating  shaft 
and  the  stationary  outer  race.  If  the  shaft  moves  in  either  direction 
the  pin  wedges  tighter.  If  movement  is  desired  in  one  direction  the 
lever  pushes  the  pin  into  the  larger  annular  space  and  allows  move¬ 
ment  only  in  this  direction.  This  type  of  a  lock  requires  that  the 
contact  surfaces  be  hardened,  otherwise  the  very  large  forces  at 
the  contact  points  will  cause  the  surfaces  to  dimple  and  the  lock 
will  be  ineffective.  The  major  shortcoming  of  this  lock  is  that  only 
two  pins  are  used.  This  produces  radial  forces  on  the  shaft  which 
must  be  balanced  by  the  bearing  surfaces.  Under  these  high  forces 
the  bearing  surfaces  wore  and  created  much  "slop"  in  the  operation 
of  the  lock.  A  series  of  wedge  pins  spaced  circumferentially 
around  the  shaft  would  take  the  load  off  the  bearings  and  probably 
work  better.  A  device  of  this  configuration  was  found  to  be  com¬ 
mercially  available.  It  is  a  Torrington  drawn  cup  roller  clutch.  A 
model  RCBt-081214  roller  clutch  was  supplied  by  the  manufacturer. 
This  unit  which  has  a  73  in-lb.  capacity  weighs  only  .  043  lbs. 
Including  two  roller  bearings,  the  over-all  dimensions  are  3/4  in. 
diameter  by  7/8  in.  long  and  it  fits  a  hardened  1/2  in.  diameter 
shaft. 


See  Appendix. 
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In  its  purchased  state  the  roller  clutch  may  not  be 
released.  It  was  therefore  modified  as  indicated  in  Figure  V-14. 

In  its  normal  operation,  the  cage  is  kept  from  rotating  with  respect 
to  the  outer  race  by  means  of  small  protrusions  on  the  cage  which 
fit  into  grooves  in  the  outer  race.  These  protrusions  on  the  cage 
were  machined  pff,  thereby  allowing  the  cage  to  be  rotated  with 
respect  to  the  outer  race.  To  disengage  the  lock  then,  the  cage  is 
rotated  clockwise,  forcing  the  pins  into  the  larger  annular  space. 
This  prevents  the  pins  from  wedging  and  allows  the  shaft  to  rotate 
in  either  direction.  Bidirectional  locking  can  be  obtained  by 
mounting  another  unit  on  the  shaft  which  locks  in  the  opposite 
direction.  The  cage  itself  could  not  be  manipulated  externally  so  a 
bushing  was  expoxied  to  the  cage  and  provided  an  external  connec¬ 
tion  for  a  solenoid  actuator.  This  lock  was  mounted  on  the  elbow 
shaft  of  the  experimental  apparatus  and  was  found  to  work  very 
well.  The  lock  has  supported  7.  5  ft-lbs  and  can  be  released  by 
turning  the  bushing  a  few  degrees.  A  small  solenoid  is  used  to 
produce  this  rotation.  This  is  possible  because,  in  use,  the  lock  is 
only  released  when  it  is  unloaded.  That  is,  the  lock  is  only  re¬ 
leased  when  the  muscle  is  fully  supporting  the  load.  The  criteria 
for  turning  the  lock  on  or  off  will  now  be  discussed. 


FIG.  3ZT-I4  RELEASE  PRINCIPLE  OF  ROLLER 

CLUTCH 
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The  easiest  way  to  control  the  lock,  in  terms  of  hardware 
at  least,  is  an  on-off  switch.  The  patient  would  then  be  required 
to  control  the  position  of  his  limb  and  the  state  of  the  lock.  This 
may  be  acceptable  if  he  were  expected  to  control  only  a  single  joint 
but  it  would  place  a  large  burden  on  the  patient  if  he  were  to  con¬ 
trol  a  multipoint  orthosis.  The  decision  was  made  therefore  to 
control  the  actuation  and  de-actuation  of  the  lock  automatically. 
Different  schemes  would  be  necessary,  of  course,  depending  upon 
whether  the  position  or  velocity  system  were  used-  A  lock  control 
system  was  first  of  all  synthesized  for  the  position  system.  This 
system  is  shown  schematically  in  Figure  V-15.  In  principle  there 
are  actually  three  control  systems  contained  within  this  system:  a 
position  control,  a  zero  torque  control,  and  a  lock  actuation  con¬ 
trol.  The  function  of  the  zero  torque  loop  is  to  reduce  the  torque 
on  the  brake  to  zerp  by  stimulating  the  muscle  sufficiently  to 
support  the  load.  The  lock  is  then  released  and  simultaneously  the 
control  mode  is  changed  to  position  control. 

More  specifically,  the  switches  (A,  B,  C)  are  controlled 

by  the  variables  e(error),  T  (torque  on  the  brake)  and  <j>  (velocity 

Jj  c 

of  the  output  angle).  Relays  (l,  2,  3,  4),  "and"  gates,  and  flip 

flops  F  ,  F  ,  are  used  to  perform  the  switching  lpgic.  The  table 
1  2 

in  Figure  V-l6  shows  the  switching  operations  at  times  tQ  to  t^ 


FIG. re- 16  STEP  RESPONSE  OF  POSITION 
CONTROLLER  UTILIZING  A  JOINT  LOCK 
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during  the  response  to  a  positive  command  signal  (  A^^).  In  the 

table  logical  "l"  implies  a  signal  of  sufficient  magnitude  to  activate 

the  relays,  or  in  the  case  of  switches,  a  switch  closure. 

At  t^  the  system  is  at  rest.  Switch  A  is  closed  making  it 

a  closed  loop  position  system.  Switch  C  is  also  closed  so  that  the 

lock  is  actuated  and  imposes  a  high  stiffness  (K)  to  the  joint.  At 

t  ,  e  >  e  j ,  and  a  stimulus  I  is  applied  t<?  the  muscle.  This 

value  of  e  also  passes  through  relay  1  then  back  to  flip  flop  F^, 

and  causes  switch  B  to  close  and  A  to  open.  The  system  is  now 

a  closed  loop  torque  system  with  zero  input.  The  stimulated  muscle 

causes  T^  to  increase  thereby  relieving  the  joint  lock  of  the  load 

torque  (T  )•  When  the  brake  torque  (T  )  is  "0U  this  signal  passes 
’  L  x5 

9 

through  3  and  is  "anded"  with  e>e1  and  <f>c  =  0  and  causes  F2  to 
switch  to  "0".  When  this  happens  (t^)  two  things  occur  simul¬ 
taneously:  switch  C  opens  releasing  the  brake  and  F^  is  switched 
causing  B  to  open  and  A  to  close.  The  system  is  once  again  a 
closed  loop  position  control  and  <j>c  increases.  The  value  of  e 
reaches  e^,  at  t^  and  the  stimulus  I  goes  to  zero.  When  e,  T^, 

and  <f>  =  0,  F_  switches  to  1  closing  C  which  actuates  the  lock 

vc  2 

(t  ).  The  system  is  now  at  rest  in  the  new  position.  In  a  similar 
manner  the  response  to  a  negative  command  signal  can  be  described. 
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The  closed  loop  position  system  as  we  have  previously 
seen  is  very  load  sensitive,  and  a  complex  compensation  scheme 
would  be  necessary  to  make  it  at  all  feasible. 

A  simpler  solution  was  therefore  sought  using  the  velocity 
type  of  controller.  The  system  to  be  discussed  was  developed 
jointly  by  the  author  and  Barry  Romich,  undergraduate  student,  for 
his  senior  lab  project.  After  many  modifications  of  the  original 
design,  the  system  shown  in  Figure  V-17  was  finally  implemented. 
This  system  also  contains  three  different  control  systems:  an  open 
loop  torque  coptrol,  velocity  control,  ?md  a  digital  control  for  the 
lock.  At  rest  the  lock  is  actuated  and  supports  the  load.  To  begin 
any  movement  the  input  signal  must  be  raised  until  the  muscle 
supports  the  load.  This  is  the  open  loop  torque  system.  As  the 
input  increases  to  support  the  load  it  triggers  the  flip  flop  which 
energizes  a  solenoid  to  release  the  lock.  Because  of  the  character¬ 
istics  of  the  lock  it  may  not  release  until  the  torque  on  the  lock  is 
zero.  The  lock  therefore  releases  only  when  the  muscle  fully 
supports  the  load.  At  this  stage  the  limb  may  be  moved  in  a  posi¬ 
tive  or  negative  direction.  The  lock  is  then  actuated  by  a  slight 
positive  velocity  (against  the  load),  or  in  the  case  of  a  movement 
in  the  negative  direction  (with  the  load),  when  the  limb  comes  to 
rest.  The  lock  always  allows  movement  in  the  positive  direction. 
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This  complete  system  was  implemented  as  shown  in  the  schematic 
diagram  of  Figure  V-18.  More  complete  circuit  diagrams  may  be 
found  in  the  Appendix,  Notice  that  the  three  types  of  velocity  con¬ 
trollers  may  be  used:  open  loop,  proportional  closed  loop,  and 
proportional  plus  integral  closed  loop. 

3.  Evaluation  of  Velocity  Controllers 

Experiments  were  conducted  using  the  above  described 
control  systems.  In  each  case  a  permanent  record  was  made  on  a 
Visicorder  recorder  of  the  command  signal  (</>r)»  the  stimulus  (i), 
the  angular  velocity  of  the  elbow  (d^)  and  the  elbow  angle  itself 
(d  ).  (See  Figure  V-19).  The  first  point  of  interest  was  the 
stability  of  the  proportional  closed  loop  velocity  system.  The  loop 
gain  was  increased  until  the  velocity  oscillated  considerably  when 
an  input  signal  was  applied.  Figure  V-20  shows  this  underdamped 
response.  Notice  that  the  frequency  of  oscillation  is  approximately 
3.  0  cps.  This  value  is  very  close  to  the  frequency  predicted  by  the 
MIMIC  simulation  (^2.  94  cps)  and  the  value  calculated  in  Section 
VB1  (20  rad/sec  =  3.  18  cps).  The  critical  value  of  loop  gain  was 
also  found  in  this  section  to  be  3.  15.  The  gains  of  the  electronic 
components  in  the  loop  are  easily  available  from  Figure  V-18.  The 
operational  amplifier  M-2  was  set  at  a  gain  of  3.  5  volts /volt;  the 
modulator  had  a  gain  of  1.  7  rpa/volt,  and  the  tachometer  had  a 


176 


2  cc 


FIG.  3C- 18  SCHEMATIC  DIAGRAM  OF  VELOCITY  CONTROLLER  UTILIZING  JOINT  LOCK 


STIMULUS  I  -  I - -  4- 


FIG.  V -20.  UNDERDAMPED  RESPONSE  OF  PROPORTIONAL 
CLOSED  LOOP  VELOCITY  CONTROL. 


FIG.  V -21.  UNDERDAMPED  RESPONSE  OF  PROPORTIONAL 
PLUS  INTEGRAL  CLOSED  LOOP  VELOCITY  CON¬ 
TROL. 
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sensitivity  of  .  025  volts /rad/sec  which  was  amplified  by  a  factor 
of  100,  The  loop  gain  was  then 


3.  5  x  1. 7  x  x  x  2.5  =  3.  15 


(V.  20) 


which  is  a  reasonable  value  for  K  as  determiped  from  the  curve 
of  Figure  III-27. 

The  same  procedure  was  repeated  for  the  proportional  plus 
integral  controller.  The  underdamped  response  is  shown  in 
Figure  V-21.  The  frequency  of  oscillation  is  2.08  cps.  The  value 
predicted  by  the  MIMIC  simulation  is  2.  27  and  the  value  calculated 
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in  Section  VB  1  is  15  rad/sec  -  2.38  cps. 

The  gains  of  the  closed  loop  systems  were  then  set  to  give  a 
slightly  oscillatory  response.  The  three  systems  were  then  sub- 
jected  to  step  changes  in  inpnt  The  average  value  of  the  out¬ 

put  velocity  (<J  )  was  then  recorded.  The  data  is  plotted  in  Figure 
V-22.  Notice  that  the  straight  line  representing  the  open  loop 
operation  does  not  go  through  the  origin.  This  is  due  to  the  biasing 
of  the  operating  point  by  the  load.  Increasing  the  load  causes  this 
line  to  shift  down  but  remain  at  the  same  inclination.  Before  a 
movement  can  be  performed  the  input  must  be  raised  to  a  value  of 
1  volt  (in  this  case)  before  the  load  torque  is  balanced.  Closing  the 
proportional  velocity  loop  has  no  effect  on  this  bias  level  since  the 
system  is  actually  running  open  loop  qntil  this  value  is  exceeded. 
The  proportional  plus  integral  controller  on  the  other  hand  does  not 
have  a  biased  operating  point.  Any  positive  nonzero  input  will 
cause  the  limb  to  move  at  a  qontrolled  velocity. 

Knowing  the  open  loop  response,  ^he  closed  loop  response 
can  be  predicted.  That  is,  using  the  known  values  of  gain  and  the 
muscle  model  the  system  equations  can  be  written.  For  example, 
the  open  loop  response  for  a  gain  of  8.  50  ma.7  volt  can  be  fitted  by 
the  straight  line, 

l  =  1. 70  I  -  1. 8 
^  c  r  r 


(V.  23) 


OUTPUT  -  (RAD/SEC) 
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FIG.3E  -  22,  INPUT-OUTPUT  RELATIONS  OF  OPEN 

AND  CLOSED  LOOP  SYSTEMS 


This  line  is  drawn  in  Figure  V-22.  It  is  seen  to  represent  the 
data  points  very  well- 
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Similarly,  for  the  proportional  plus  integral  controller; 
assuming  the  muscle  dynamics  to  be  the  same,  the  steady  state 
response  is  given  by 


'f 


lim  _ 

3-^0  , 

s(. 


2(1.  7s  +  10.2)(.  lie"’025) 


05s  +  1)2(.  15s  +  l)  +2  (1.  7(3+10.  2)(.  116‘  °2s)2.  5  J 


(V.  29) 


-  (2)(10.  2)(.  11)  rad /sec 

(2)(10. 2)(.  11)(2.  5)  *  volt 


(V.  30) 


The  slope  of  the  line  through  the  experimentally  obtained  points  of 
Figure  V-22  is  0.  55  rad/sec/volt.  A  dotted  line  of  slope  0.4 
rad/sec/volt  is  also  drawn  on  the  figure  for  comparison. 

The  preceding  discussion  has  shown  that,  from  the  point  of 
view  of  control  theory,  the  controllers  perform  pretty  much  as 
expected.  What  needs  to  be  evaluated  next  is  whether  the  lock  and 
its  associated  logic  circuitry  works  and  if  so  can  a  human  use  the 
controller.  Practice  sessions  were  held  to  see  how  well  the  posi¬ 
tion  of  the  arm  could  be  positioned  using  each  of  the  controllers. 
The  lock  actuation  control  was  found  to  work  quite  well  except 
for  lock  actuation  when  the  limb  has  a  negative  velocity.  In  this 
condition  the  lock  should  then  actuate  when  the  velocity  of  the  limb 
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approaches  zero.  This  mechanism  was  found  to  be  unreliable  and 
would  cause  the  lock  to  actuate  prematurely.  In  any  case  the  lock 
may  still  be  actuated  by  a  positive  velocity  of  the  limb.  That  is, 
movements  in  the  negative  direction  must  be  followed  by  a  slight 
positive  velocity  to  actuate  the  lock.  With  the  lock  control  system 
performing  satisfactorily  then  the  velocity  controllers  were  ap¬ 
praised.  The  desired  response  independent  of  the  type  of  con¬ 
troller  was  the  same.  Figures  V-23  and  24  show  successful 
attempts  to  position  the  limb  in  both  the  positive  and  negative 
directions.  Notice  in  Figure  V-23  that  the  input  ({_;)  is  brought  up 
to  a  value,  the  arm  flexes  at  a  constant  rate,  and  when  the  desired 
position  is  reached  the  input  is  reduced  rapidly  aud  the  limb  stays 
locked  in  this  position.  In  Figure  V-24  the  arm  extends  to  a  posi¬ 
tion.  This  is  brought  about  by  first  of  all  increasing  the  input  until 
the  muscle  overcomes  the  load.  At  this  point  the  lock  releases  and 
the  input  signal  may  decrease  causing  the  limb  to  extend  with  the 
load.  The  desired  angle  is  overshot  a  slight  amount  then  the  input 
signal  is  once  again  increased  to  overcome  the  load  and  cause  the 
lock  to  actuate,  then  decreased  to  zero.  The  relative  difficulty 
with  which  these  operations  can  be  performed  with  each  of  the 
three  velocity  controllers  is  discussed  below. 


STIMULUS 


FIG.  V -23.  CONTROLLED  FLEXION  USING  PROPORTIONAL 
CLOSED  LOOP  SYSTEM. 


STIMULUS 


FIG.  V-24.  CONTROLLED  EXTENSION  USING  PROPORTIONAL 
CLOSED  LOOP  SYSTEM. 
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Open  loop  velocity  controller.  -  This  system,  as  we  have 
seen,  ,  is  very  load  sensitive.  This,  however,  may  not  be  as  large 
a  drawback  as  thought  earlier  since  the  human  controlling  the  limb 
is  quite  adaptive.  If  a  patient  were  to  control  a  single  muscle  he 
may  be  able  to  do  this  with  no  feedback  fairly  well.  In  the  experi¬ 
ments  the  open  loop  gain  was  fairly  high  which  means  that  inputs 
that  were  slightly  above  the  equilibrium  point  caused  the  limb  to 
move  fairly  rapidly.  Even  so,  a  fair  degree  of  control  was 
achieved. 

Closed  loop  proportional  velocity  controller.  -  Feeding 
back  the  velocity  does  not  affect  the  amount  of  input  signal  needed 
to  reach  the  equilibrium  point  before  movement  takes  place.  This 
means  that,  as  in  the  open  loop  case,  the  human  must  "hunt"  for 
this  value  of  input  before  he  can  perform  a  movement.  Above  this 
value  the  velocity  is  more  easily  controlled  than  the  open  loop  since 
there  is  some  regulation.  Again  the  adaptive  behavior  of  the  human 
must  be  exploited  for  this  system  to  perform  properly. 

Closed  loop  proportional  plus  integral  velocity  controller,  - 
This  controller  proved  to  be  superior  to  the  other  two  in  that  there 
was  no  "dead  zone"  for  the  input  signal.  There  is  no  hunting  for 
the  operating  point  and  any  positive  signal  causes  the  muscle  to 
contract.  There  is  a  lag,  however,  due  to  the  integrator  in  the 
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forward  loop.  After  some  training,  however,  the  adaptive  human 
is  able  to  use  this  controller  as  well  or  better  than  the  other  two. 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusion  is  drawn  that  controlled  stimulation  of 
muscle  is  indeed  feasible  for  the  actuation  of  limb  movement.  As 
mentioned  previously,  however,  all  experiments  were  performed 
on  normal  relaxed  individuals,  simulating  muscle  with  an  upper 
motor  neuron  lesion.  As  far  as  the  stimulated  unparalyzed  biceps 
muscle  is  concerned,  it  produces  enough  torque  (about  4  ft-lbs)  to 
power  the  limb  and  a  modest  load  against  the  force  of  gravity. 
Further,  these  movements  can  be  performed  at  speeds  within  the 
range  for  normal  movements  (up  to  6  rad/sec).  The  controllability 
of  these  movements  using  either  of  the  velocity  controllers  dis¬ 
cussed  and  a  joint  lock  is  also  quite  good.  The  proportional  plus 
integral  controller ,  though,  seems  to  be  superior. 

It  seems  that  there  is  no  functional  difference,  in  terms  of 
stimulated  response,  between  normal  relaxed  muscle  and  muscle 
with  an  upper  motor  neuron  lesion  (UMNL).  The  dynamic  modeling 
and  control  schemes  therefore  should  be  directly  applicable  to 
non-atrophied  muscle  with  an  UMNL.  The  parameter  values  of  the 
muscle  model  will  change  from  one  muscle  to  the  next,  but  the 
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form  of  the  model  should  be  directly  applicable  to  all  of  the 
muscles  of  the  body.  It  is  very  unlikely,  however,  that  a  paralyzed 
person  can  be  found  whose  paralyzed  muscles  all  have  an  UMNL. 
Paralyzed  persons  generally  have  both  UMNL  and  LMNL's.  The 
electrical  excitability  of  muscle  with  a  LMNL  has  not  been  dealt 
with  in  this  thesis  other  than  to  indicate  the  form  of  its  strength 
duration  curve.  If  we  consider  the  optimal  pulse  width  to  be  the 
chronaxie  value  (shown  by  Vodovnik  in  ref.  67),  then  the  optimal 
pulse  width  for  muscle  with  a  LMNL  is  about  2  msec.  This  means 
that  the  duty  cycle  of  the  DC  pulsed  stimulus  must  be  increased  by 
a  factor  of  10  (0.2  msec  =  chronaxie  for  UMNL).  This  increase  in 
the  energy  level  of  the  stimulus  will,  of  course,  be  more  painful 
as  indicated  in  Figure  111-8.  Perhaps  a  different  waveform  would 
be  a  better  stimulus. 

In  any  event,  any  comprehensive  orthotic  aid  powered  by 
stimulated  muscle  would  involve  two  different  types  of  electrical 
stimuli  to  actuate  the  differently  involved  muscles.  It  is  therefore 
recommended  that  the  mechanical  response  of  electrically  stimu- 
lated  muscle  with  a  LMNL  also  be  modeled  and  evaluated  as  to  its 
feasibility  as  a  controlled  actuator. 

The  solutions  to  controlling  a  stimulated  limb  described  in 
this  thesis,  though  workable  and  quite  effective,  require. considerable 
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instrumentation  both  on  the  orthosis  and  off.  Since  this  system  was 

built  and  tested,  however,  new  insights  into  the  problem  have 

emerged  and  it  now  seems  possible  that  a  much  simpler  controller 

could  be  built  which  may  work  as  well  or  better.  This  system  is 

the  closed  loop  torque  control  shown  in  Figure  VI- 1  .  The  torque 

under  control  (T  )  is  the  resultant  of  the  muscle  torque  minus  the 
c 

load  torque.  This  is  the  actual  torque  available  to  move  the  limb 
which  has  a  fixed  gain  of  about  2  rad/sec/ft-lb,  depending  upon  the 
viscous  damping.  This  torque  can  be  transduced  by  placing  strain 
gages  on  the  splint  structure.  Controlling  this  torque,  therefore, 
controls  the  output  velocity  of  the  limb.  The  control  can  be 
effected  by  varying  the  amount  of  torque  feedback  (H).  At  rest 
H  =  0  and  the  stimulus  is  off.  To  begin  a  movement  H  is  set  to 
1. 0  and  the  muscle  produces  a  torque  equal  to  the  load  (T^  =  0). 

This  eliminates  "hunting"  for  the  zero  torque.  The  lock  is  then 
released  using  the  same  switching  logic  as  the  velocity  controllers. 
Varying  H  about  the  1.  0  value  then  controls  the  velocity  of  the 
limb . 

Implementation  of  this  scheme  requires  no  bulky  potentio¬ 
meters  or  tachometer  to  be  mounted  on  the  splint  structure.  All 
that  is  required  is  a  joint  lock  at  the  bearing  and  the  application  of 
strain  gages.  The  electronic  hardware  necessary  is  a  high  gain 
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amplifier  with  variable  gain,  a  modulated  stimulator  and  a  compen¬ 
sation  scheme  to  allow  the  necessarily  high  loop  gain.  It  is  recom¬ 
mended,  therefore,  that  this  controller  also  be  implemented  and 
evaluated  along  with  the  velocity  controllers.  A  lightweight  ortho¬ 
sis  should  be  developed  which  includes  the  joint  lock  and  control 
schemes.  Then  by  a  carefully  conducted  evaluation  with  patients 
this  device  could  be  directly  compared  with  existing  orthoses 
powered  by  motors  and  cross-weave  actuators . 

It  may  be  found  that  tracking  a  motor  point  with  a  surface 
electrode  is  not  a  practical  solution.  In  this  case,  then  implant 
stimulation,  where  the  stimulating  electrode  is  placed  on  the 
muscle  or  nerve,  may  be  the  answer.  Work  is  being  carried  out  in 
this  area  by  members  of  the  Cybernetic  Systems  for  the  Disabled 
group.  Dr.  Vodovnik  and  graduate  students,  R.  Lorig  and 
H.  Peckham,  are  investigating  the  use  of  radio  frequency  telemetry 
as  a  control  signal  for  an  implanted  stimulator. 
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APPENDIX 

GRAPHICAL  DECONVOLUTION 


Problem  statement  -  Given  the  transfer  function  and  output  of  a 
blackbox,  determine  the  input. 

Method  of  Solution  -  The  input  will  be  approximated  by  a  sequence 
of  impulses  of  various  magnitude  such  that  the  sequence  of  impulse 
responses  at  the  output  add  to  give  the  measured  output  curve. 

Once  the  input  impulses  are  determined,  a  smooth  input  function 
can  be  synthesized. 

Procedure  - 

1.  On  transparent  graph  paper,  construct  a  set  of 
impulse  response  curves  for  the  black  box  with 
different  magnitudes .  (Figure  A- 1 . ) 

2.  Plot  the  output  curve  on  similar  graph  paper 
using  the  same  time  scale.  (Figure  A-2.) 

3.  Now  place  the  graph  paper  with  the  output  curve 
over  the;  impulse  response  curves  so  that  the  time 
axes  coincide.  Slide  the  output  curve  to  the  right 
or  left  until  the  output  curve  most  closely  matches 
9.n  impulse  response  curve. 

4.  Now  subtract  the  impulse  response  curve  from  the 
output  curve.  Record  the  magnitude  of  the  impulse 
on  the  output  plot  at  time  t=  0  of  the  impulse  response 
curve.  (Figure  A-3.) 


Figure  A-l.  IMPULSE  RESPONSE  CURVES  OF  BLACK  BOX 


OUTPUT  CURVE 


I 
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5.  Repeat  this  procedure  on  the  reduced  output  curve 
until  this  curve  has  a  near  zero  value  for  all  time. 

The  resulting  sequence  of  impulses  then  repre septs 
an  approximation  to  the  input  curve.  If  these 
impulses  are  evenly  spaced,  then  a  smooth  curve 
may  be  drawn  through  the  peaks  of  the  impulses. 

If  the  impulses  are  not  evenly  spaced,  an  equivalent 
set  of  evenly  spaced  impulses  may  be  obtained. 

Integrate  the  pulses,  then  fit  the  resulting  stepwise 
curve  with  a  smooth  line.  Now  divide  the  time  axis 
into  equal  increments  to  obtain  the  set  of  evenly 
spaced  impulses. 

The  procedure  is  quite  rapid  and  allows  a  person  to  guess  at 
a  solution  and  try  it  immediately. 
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Figure  A -4.  TORQUE  TRANSDUCER 
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Part  Weight  (lbs)  Inertia  (Ft  lb  sec  ) 


1 

.  88 

. 02090 

2 

.44 

. 00096 

3 

.  67 

. 00087 

4 

.46 

. 00012 

5 

3.  52 

.00009 

6 

1. 00 

. 00013 

ITotal 

. 02307 

Figure  A -5.  INERTIA  OF  ARM  SUPPORT 

STRUCTURE 
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Assuming  a  coefficient  of  friction  y  =  0.2, 
F1  “  MP2-:-0.2P2 

For  a  torque  of  T  =  5  ft-lbs,  and  a  shaft 
radius  r  =  .25  inch, 

F1  -  I  “  tUt  -  240  lbs 


F2  =  1200  lbs 


In  the  multi-pin  lock,  this  force  is 
distributed  among  the  pins  and  reduces  the 
dimpling  effect  of  the  pins  on  the  shaft 
and  outer  race. 


Figure  A-6.  FORCES  IN  ORTHOTIC  LOCK 
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